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SYl'TOPSIS 


ANALYSIS OP RSILPORCEI) 
CUj‘ CRETE CYLILOra-CAL SHELL BDOES 


(a thesiD cuLmit bed in partial fulj.il 'cent of tne requirement 
for the decree of Doctor of Philosophy hy K. Surayya to the 
Departi'ienb ol Civil Dnyineerindj Indian Institute of Techno- 
logy, Kfi.-nurj June 1978) 

A structure may he analysed either hy elastic or limit 
analysis. The present study is focussed on the limit analysis 
of reinforced concrete shell roofs, Luait analysis is based 
on two fundamental theorems, viz.,(i) static (or lov/er hound) 
theorem oiid (ii) kmematic (or upper hound) theoreui. These 
tv/o thcoreras provide lovirer and upper hounds to the actual 
collapse loud- 

I'i.If. Eialkow proposed both lower and upper hound 
solutions for simply supported cylindrical shell roofs. Rigid 
plastic Uxialysis with linear approximation to the Hises and 
Tresco yield criteria is adopted, II, A. Snoeh and ¥. C. 
Schnohricu developed elast ic— plastic solution for cylindrical 
shell roofs using the Mises yield criterion. A.L,L, Baker 
considorod circular reinforced concrete cyl ■Uidr ical shell as 
a craclcoj hoam of curved cross sectioxi. M, Janas atid A.E. 
Rshanits^n jproposed upper hound solutions independently, for 
simnly supported reinforced ccaapxete cylindrical shell roofs. 
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C .S . Lin 'ni'j A.C. Scordelis applied finite element bechniq_ue 
to analycG reinforced concrete shells of general form, 

llif present study is concerned with the development 
of lower Loi'nd solutions for simply supported reinforced 
concrete circilar cylindrical shell roofs subjected to uni- 
form gr? viby loading. The longitudinal edges of blie shells 
considered liercin are either free or with edge heaiiic, 

tv luiiMUC yield condition for a reinforced concrete 

cylindnco] shciH element is not available in the literature, 

Laced on coi.e nccumptions b\/o approximate yield criteria are 

pro[)osod, ^1.. .,(i) - yield criterion end (ii) - ^4, 

yield cribo-’ion. The H -yield criterion states tliab the 

c 

choll cL.r"'ib yields if the longitudinal compression, 
roaches leiuc, U , v/hich is equal to the strength of tho 
shell in compression per unit length. The - h yield 
criterion r; bates tnat tho shell elenignb yields duo to the 
intoracbioii 'n bending moment, M^, and inplane direct force, 
ill b'lc croneverse direction. In this study yielding of 
the shell olcmcnt is defined if at least one of the two 
yield co.'io'ons is satisfied. To simplify the analysis, the 
M(j, - ^/idd criterion v/hich is nonlinear is aporo^cimated as 
pioce-v/ipo linear- 

In general a cylindrical shell element is subjected 
tx) inplanc direct forces end , inplane shears and 



XV13- 


■trc'-nrjv'Tse shears and Q^, hending manents II^ and , 
and twisuing moments M and M * But for thin shells they 
are reduced to Q^, Q^, and In this 

study, tuc 0 tress resultants CL, M and 11 are neglected. 

X. }C X 9 

Since these stress resultants cannot be neglected for short 
shells, b’-'csc solutions are not applicable to this class 
oC shells, lie equilibrium equations are reduced to four, 
containin , :.3--ve unkno^m stress resultants, 

llsni,; static theorem lower bound solutions are deve_ 
Jopod. liiico i,hc number of stress resultants are nore than 
the equ Llibr lU ' 1 equations by one, N distribution is assumed 
in all t.ic r.'"lations. Its magnitude is determined from the 
static Oi'u l 1 ibrium conditions of the shell, Knoi/mg 11^, 
other strcca resultants are solved using equilibrium equations 
and boundory conditions for the stress resultants. Three 
solutions are developed vfith three different distributions 
for B . In the first solution, N is taloen as uniform both 
in tension and coiopression zones across the section. In the 
second solubion, iLr vax^es linearly in tension zone aixd uni- 
form in (;oumrG&oion zone. In the third solution, parabolic 
distribution as tsUcen in the compression zone keen mg linear 
variation an tension zone. In the solution three, ^ the 
angle defining the position of maximum longitudinal compre- 
ssion IS kOot as a variable to obtain the optimum lower bounds. 
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For shells itn edge beams, the longitudinal force developed 
in the edge beam is taVen as trapezoidal. 

Shells '/L'cli different geometric parameters are analy- 
sed by boon clastic (aSCE, kanual No. 3i) and by tne proposed 
lov'/er bound solutions. Distribution of stress resultants at 
critical sections are presented for comparative st ■'dj?'. Two 
design examples are also presented by elastic analysis for 
service load and by the proposed limit analysis for a load 
factor of The economy of the materials is discussed. 

A nui-nrical technique ic also developed. In this, 
tne problem is formulated as a linear programming problem and 
tbe revised simplex method is used to solve it. The shell is 
considered as a grid of discrete points. Ninite difference 
scheme is used for equilihr lum equations. The load is tahen 
as an objective function to be optimized and the stress resul- 
tants at the grid points are taloen as variables- The con- 
straints are formulated in terms of these variables at each 
grid point to satisfy the yield condition. The load is opti- 
mised usinu 'die simplex method and bhe corresponding variables 
are evalUM-tcd. The distribution of stress resultants at 
critical fjccbionsis also presented. 

The UiGsis concludes with the discussion of the 


results dtoi^.iaed* 



CHAPl'^R 1 


Il'JTROnJCTIOJ’4 

1*1 General 

The \, ord ^ shell’ primarily means any curved shape/ 
surface endosi/jg/cover ing, certain volurae/area. It nas 
different connotations for persons in varied enviromicnt . 

The shells \/o find around in abundance are either God's 
creation or uaii—made. The hard outside covering of a seed, 
a fruit or an '’nimalj the ha.rd calcareous envelope ol a 
bird's egg [all m the former category. An explosive pro- 
;]ectLle uc'd mi a cannon or a mortar may be categorised in 
the latter f riety. If we stuny the nature in deu oh it will 
at once bo realised that microscopic cello in tlio living 
beings, root sections of birds' feathers aiid the hoinboos 
used in building construction from the times innenov ..ii arc 
also shell shaped . The egg &..clls range widely ii. sire and 
shape . 

Linn by observing narure critically, realised tact a 
tgivcn quantum of material ii‘ lai3 out judiciously con possess 
tremendous looui carrying capabilities. The ob-vious shape is 
tuat OJ a shell. Hov/ virtually every item of the most 
modern industrial equipment, he il aero-space, bulletin,,, 
nuclear, marine or petro-chemical, has shell as a structural 
comp onent « 



all structures are invariably three— 
dimoxisional, lov purposes of analysis they are charactorised 
as one-diiaenn j-oxial, tv/ o- dimensional or three-dimensional, 
depending luo.l Mie relative ratios of their principal dimen- 
sions. In biiat: sexise, shells are tv/o-dimensional (thickness 
IS small in conparicon '‘/ith the ocher two principal dimensions) 
and are represented by the middle surface. Thus an arbitrary 
point of the middle surface of a shell may be defined by 
specifying fcuo scalar components. The complete geometry of 
a shell IS defined by describing the middle surface and the 
thiclcnoss . 

As defined earlier, a shell is a curved surface. It 
is called thixi if the maximum value of the ratio of the 
thickness tp the radius of curvature is very small in 
comparison with unity. The shells may be curved in one or 
two principal directions and are termed singly or doubly 
curved resuectively, A cylindrical shell falls in the former 
category while a spnerical one falls in the latter. 

In tills thesis a cylindrical shell is defined as 
a curved slab (Gin,_^ly curved shell) cut out from a cylinder. 
The slab is thus bounded by tv/o Gtrai,_iht longitudinal edges 
parallel to t/io axis of the cylinder and by tv/o curved 
transverse odges in planes perpendicular to the axis. The 
longitudinal edges are either free or stiffened with beams. 
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ShellG nay “be made of difxereni; materials depending 
uoon tno luncbJonal, strength and weight requirements. In 
building inCisfcry most of tiie shells are made of reinforced 
concrete, hcmiorcod concrete shell structures possess many 
of the qualities of an ideal structure. They require a 
minimum, amount oi material and permit coverage of large areas 
with out intermediate supports. Their form offer virtually 
unlimited possibilities for aesthetic and architectural 
considerations* Cylindrical shells are also popularly 
adopted for their simplicity in form work. 

Rigor olig and approximate analytical methods of 
analysis and dorjign are readily available in literature. 
However these ore mostly confined to elastic theory. It has 
long been recognised that reinforced concrete shells sire 
usually too heavily loaded to remain in the range of validity 
of linOcir clastic theory. In addition this method of analysis 
has one more limitation that the behaviour of a structure 
can he studied within a sp® cified limit of v/orlcing stress. 

The factor of safety here is v/ibh respect to tlie ultimate 
stresses oC t]ic materials but not ’./ith respect to ultimate 
or collapse load of a structure. In practice elastically 
designed structures ar© conservative resulting in wastage 
of precious nutorials. Ilencc rational methods of analysis 
have come into vogue talcing the collapse loads into consi- 
deration. Tnese are known as plastic limit or simply limit 
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analysis rac L-.'oris, v/hich take into account the true factor of 
safety against collopse. fhe ratio of collapse load to 
working loaO ns known as load factor. Every method has its 
ovm limitabionr; . The limit analysis methods are no exception. 

The "behaviour 0i the structure at v/orking load cannot be 
predicted by ohcce methods. The present studi is confined 
to the limit anedysis of reinforced concrete cylindrical 
shell roofs. In the follov/ing paragraph a comprehensive 
literature s.urvoy is presented. 

1,P Select IV'- literature Survey 

The nro^rious investigations on the limit analysis of 
shells have restricted mainly to axis ymaetric snolls 

(l,K,3,4)^ c;:c<-ot for a fo’ • on cylindrical shell roofs, 

Pialkov'/ (5) pcoiosed both lower and upper bound solutions 
for simply suuoorted long cylin'lrical shell roofs with free 
longitudinal edges under radiol loading. An insight into 
the limit analysis of cylindrical shell roofs is given. 

Tne matej’iel is idealised as a rigid plastic material obeying 
the Tresca ox linearized !ii&cs 5 /ield criterion. Slioeb and 
Schnohrich (G) presented a solution for elastic-plastic 
onalysis of cylindrical shell rools './ibh n discrete I'lithcma- 
tical model iic wy the PranJti-Ileuss incremental tlioory of 
plasticity logcther v'lth Iiises yield condition. c. Iang(7) 

* Furabci s in ihe brackets indicate tae reference number. 
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obtciined els o oic-plastic solution for the limit analysis of 
cylindrical shells using finite element technique. The 
material obeys 'che Mises yield condition. 

Thus tiic limit analysis is mainly confined to materials 
obeying thr ilices and Tresca yield criteria. 

Very ,.o\r upper bound solutions are available for 
rotationaliy syimictric reinforced concrete shells like spheri- 
cal domes ajid conical shells (8,9). Simple failure mechanisms 
were assumed Piid the collapse loads v/ore arrived ot by equat- 
ing the external v/ork done by the load during deformation 
to the jiiternnl energy of dissipation. Lov.er and upper 
bound solutions ere available for axis ymie trie shells like 
v/atei tatilcs (lO)* Balcer determined the ultimdbe load 

oC re mi orcecl concrete cylindrical shell roofs. 't''e shell 
\/as considorc'd as a crachcd beam oJ' thm curvjd cros;j section 
in the lont^iiudinal direction srk as an aicn in oiic transverse 
direction. This analysis was not based on toe fun' i-ntals 
of limit aiio.lysis but on the orinciplo of evac' cd r^ uiforced 
concrete section, Rziianits 3 ai ( 13) devoloocd a simplified 
kinematic method for shallow Siielis. Savveauk (l4) presented 
upper bound solutions for reinforced concrete c./lindrical 
shells resting on rigid supports with lifferent edgej condi- 
tions. These \'/ere based on colla.pne mechanisms suggested 
from experincntc, Janas ( 15) prooosed Ixmit analysis of 
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nonsymmetricc'a chells by generalized yield line me-tliod. 
Rzhanitsyn (i6) introduced a colution to deterrame the 
bearing capo city of reinforced concrete cylindrical shell 
roof supported by end diophragms . The extremum principle of 
the kmomotic method of limit equilibrium v\ras used as the 
basis for this analysis. Lin and Scordelis (l7,l8) apolied 
the Iimto eleiaciib technique to analyse reinforced concrete 
shells of gencra.l form. The analysis includes the load- 
deformation response and craclo propagation and determines the 
internal concrete jjnd steel stresses. 

BoslcIgs, some experimental results of reinforced 
concrete shell roofs are also available in the literature. 
Bouma or al. (l9) conducted a senes of tests on cyJinJrical 
shells with edge beams to demonstrate the ultimate load 
behaviour uf shells- A. scries of eleven shcils v/ith the 
same cross section but of different spans were tested. The 
load factor rouged from 2,75 to 4.50* Enami (2o) conducted 
a number of experiments on reinforced concrete folded plate 
structures end based on the failure patterns presented a 
series of me oh. 'iiisms . Harris ond 7hite (2l) conducted 
experiments on small scale models to investigate clastic and 
inelastic behaviour of reinforced concrete shell roof 
structures , 

from tile forgoing it is obvious that most of the 
structures are confined to arriving at upper bound solutions 
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which are inherently unsafe in the absence of a correct 
collapse uic cLianism, It is imperative that methods of analysis 
based on the lov/er bound theorem are more appropriate* Hence 
this study IS confined to lower bound methods of analysis. 

1.3 Object < nd Scope of Present Y/ork 

A Cew uoper bound solutions are available to determine 
the ultimate load capacity of cylindrical shell roofs. The 
major objective of the present worlo is to develop solutions 
for sim,jly suooorbed reinforced concrete cylindrical shell 
roofs based on 'static theorem’. The cylindrical shell roofs 
analjrseu in this study are classified into tvi/O groups, viz., 

(i) cylinflncal shell roofs free longitudinal edges and 

(ii) cylindric.'l shell roofs with edge beams alon,^ the 
1 ong 1 tud inal c dge c . 

PrODCiitly no unirpe j^ield condition is availa,blc in 
the libemuro for reinforced concrete cyliirl-rical shell 
element. Jxta this in view t* o independent yield conditions 
are proposed, viz.,^) '1. - yield criterion’ - yielding due 
to manimum longitudinal compression and (ii) ’M^ - il^ yield 
criterion’ yielding in the transverse direction due to the 
interaction ol the transverse moment and the inplane 
direct force I’T^ . 

Three lower hound solutions are presented for 
cylindrical shell roofs witu free longitudinal edges. Ihe 
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variation ol blie ultimate load in nondimensional form to the 
ratio of span L to the radius R is indicated by way of graphs. 
Tv7o shells are analysed with different geometric parameters 
indicating the distribution of stress resultants at collapse. 

Two shells w'ibh the same geometric parameters are designed 
for a given load factor by the proposed lower bound solutions. 
Elastic and plastic design methods are compared with material 
e conomy in view . 

for cylindrical shells v/ith edge beams the follov/ing 
aspects arc studied. Three lo^iier bound solutions are developed. 
Two shells oj. different geometric parameters are analysed 
to indicate the distribution of stress resultants at collapse. 
Two shells wLbn same geometric parameters are designed for 
a given load factor by the proposed lower bound solutions. 

The economy of the materials by both elastic and plastic 
methods is discussed. 

The ultimate load calcvilatcd from the present lov^er 
bound solutions is compared .nth the ultimate load determined 
from the upper bound solution given in the literature (22) . 

It IS found that both values are tuo same if the siiell 
fails as per the H — yield critecion. Henc^ the solution is 
exact vaien tiic shell obeys tiie h -yield criterion. 

A numerical techni(i_i'e is also developed to ariive 
at the lovfex bound solutions. It is formulated as a linear 
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progranmun^^ (LP) protDlen and revised simplex metnod is used 
to solve it, Ine shell is considered as a grid of discrete 
points. Pinite difference scheme is used to replace diffe- 
rential equations of equilibrium by algebraic equations. The 
load is taloen as an ob3ective function and the stress resul- 
tants are taicen as variables . The constraints are formulated 
in terms of these variables at each grid point so that the 
yield condition is not violated. The load is optimized using 
the simplex me chod and the corresponding variables are 
evaluated. The distribution of stress resultants is presented* 

The thesis concludes with the discussion of the results 


obtained 



aiKBTER 2 


BAWROUm THEORY 

2,1 General 

Tne nurnooe of structural analysis is to study or 
predict the recuonse or behaviour of a loaded structure. It 
IS based on .Cour requirements vis., the equilibrium, the 
condition of compatibility, the stress-strain relationship 
of the material used and the known boundary conditions. The 
equilibrium c.iiuitions specify the relationship between the 
internal streen resultants and the external loads applied. 

The relation betv/een displacements and strains are termed 
compatibility equations. The boundary conditions may be 
static (bending, moment, shear force etc.) or kinematic 
(deflection, rotation etc.) constraints on the structure and 
are usually prescribed over certain domain. The stress-strain 
relations are determined by simple tests. If the foregoing 
four conditions are satisfied by a solution at every point 
in a body, that solution is said to be exact and yields 
correct stresses, strains and displacements in the 
body. Sometimes it is not possible to satisfy all the four 
conditions, one or more of then, may be relaxed. A solution 
satisfying loss than these four conditions is approximate. 
Depending on tnc ronge of mechanical behaviour, analysis may 
be broadly classified as elastic or plastic. 
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If tile I’clati onsliip beti^/een the stresses and strains 
IS a veined (unique) and does not depend on the stress 

history (no effocts of yiscosity and hystersis) then the 
material is tf-naed elastic, '/hen the foregoing relation is 
linear then Lje uiateiial is linearly elastic (The material is 
said to obey generalised liooke’s law). Elastic anolysis 
predicts the behaviour at working loads while plastic analysis 
does at failuie . The stress-s Irain relation for a rigid- 
perfecLly plcntic material is termed flav r'ule . Because of 
the non~uni(iUC nature of the stress-strain relation for a 
plastic malo.iol» complete solutions are difficult to get 
for complen Ooructures. A-s a result approximate solutions 
known as bounds o/i the correct collapse loads are usually 
found in literature. 

In this chapter an attemot has been made to present the 
equilibrium equations of a cylindrical shell element , limit 
analysis theorems and various yield conditions. 

2,2 Equilibrium Equations 

A cylindrical shell may be defined as a curved slab 
v/hich has boo a cut out from a full cylinder. The slob is 
bounded by uio straight ’ion;',itudinal edges’ parallel to the 
axis oC the c/lmder and by two curved ’ tronoverse edges’ in 
planes perpeni.!-! culac to the axis. The slab is therefore, 
curved in oii'.y one direction. Vhen the curvature is constant, 
the cylindctca.1 snell is said to be circular. 
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The equilibrium equations are derived based on the 
followin:=, assumptions : i) The shell is thin, ii) The deflec- 
tions are snail at the instant of collapse, '.7ith the assump- 
tion (ii), the equilibrium equations are derived based on the 
undeformed feoiaetry of the shell element. 


T'lq. 2. la shOv\^D a typicol shell, its geometry and also 

the positive directions of the coordinate axes. Orifin is 

taicen at the a,pex of the shell directrix at left hand support. 

The coordinate x is measured along the crown generator, the 

coordiiiate y along the directrix, the coordinate 2 along the 

inv/ard normal. The forces and moments acting on a shell per 

unit length are usually known as the stress resultants. In 

general, a cylindrical shell element is under the action of 

inplane direct forces and , inplane shears and 

transverse shears t and , bending moments and il, and 

twisting moments M ^ and tl , The stress resultants and 

X(p <pX 

sign convention adopted are snovm in Pigs. 2.1b and c. p^, p^ 
and p^ denote the components of external load applied per 
unit area in y,y and z directions respectively. The equili- 
brium equations for the shell element are given as (23) 
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.., (2.3) 
. .. (2.4) 
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^ R = 0, 


... (2.5) 


where R j-s the radius of curvature . Ror thin shells it 
can he readily sJiown that H = N. and Li ^ . In all, 

there are five equations of equilihrium and e ight unlcnown 
stress resul hints . Thus equilibrium equations alone are 
not sufficicuiG for solving stress resultants. To rcciuce the 
complexity of Lne problem, certain simolif ications ore made 
based on tho structural action. 


In the literature some numerical studies as well as 
experimental results show that the state of stress m a shell 
depends essentially on the ratio of the span I to the chord 
v/idth 1 (24) . Depending on the ^ ratio it is possible to neglect 
certain stj .or, resultants. Accordingly cylindrical shells are 
classified es short (g- < l) , medium ( i < g- < 3) and long 
(| > 3) . 

Tho s ioit cylindi’ical shells arc charactcrisod by the 
appearance o oil stress resultants that occur m a shell. In 
the case oi ue^'ium length cylindrical shells, theoretical and 
exper ineiitni studies have shown that the scate of stress 
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IS influenced priinarily \i-y the transverse lending monent , 

while the elfecto of twistin^ moments I.I 11 , , the bending 

moment M and shear force Q_^ are negligible , Longitudinal 
X x 

shells can be treated as thin v/alled rods, based on the 
hypothesis tliot the cross section is not deformed- 


The present study is mainly concerned v^/ith medium 
cylindrical shells. It can be applied to long cylindrical 
shells as \/ell. In this thesis no distinction is made 
between medium and long cylindrical shells. 


Neglc cling the stress resultciits LI , Q and li ,, the 

x’ ^X X^)’ 

equilibrium equations (2.l) through (2.5) reduce to 


Oil 


8l! 


(J»x 


54i 


+ ^ = 0 , 


( 2 . 6 ) 


+ II ~ Q A + P K h ^ 0 , 

6<J) dx ^ <> 

oQx 

+ Pz ^ = 0 » 

1 

- E6f -*• 


... (2,7) 
. .. (2 . 8 ) 

. (2.9) 


The ciiell IS subjected to gr. vJty looding only. If p 
IS the inter. nty of loading, them its components are ,,iven by 


•a. = 0 , 

.A, 

D = p sin 6 
- 0 

p^ = p cos <> . 

Zi 




« • • 


( 2 . 10 ) 
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DifXerenbioti-i^ Lqs. (2.6) >.md (2,7) with respect to x and <P 
respectively rnd using the relations ( 2 . 10 ) yields 


5 h . d^i! 


T * dST = ° ’ 


... ( 2 . 11 ) 


I'h , 6Q, 

■w + dSdV " r/' + p a cos ♦ = 0 


... (2.12) 


Eq. (2.12) 1 C simplified ucing Eqs . (2.8) and (2.li) as 


2 sV 

Nx + — ^ R — + 2 p R cos (j = 0 • 
3 dx 


... (2.13) 


Eliminatini 


g Q from Eq. (2,9) and substituting in Eo. (2.8) 


for Q ^ one goto 


1 ^‘’^0 

4 . ^ ^ + p R cos <!) = 0 . 

^ 0 ^ 


... ( 2 . 14 ) 


Thus the couilihriums equations (2,6) through ( 2 . 9 ) 
are reduced to Coll owing equations : 


dl( Sh 

R ^ 4 . -_J 

ox 5 


. (2.15) 


6'^r, p 3 IT 

IT +■ . — ■pi _ R . — . 1 11 ^ .-. 4 - 2 p R cos 3 
3 3" 3x 

1 ^^’^3 

if + ^ -~-p— +■ p E cos 3=0, 

- r6^ + a* E= 0- 


. ( 2 . 16 ) 


. (2.17) 


( 2 . 18 ) 
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The modified form of eq_uilihrium equations (2,15) 
through (2,i8) are used in the analysis. 

2,3 Basic Concepts of Plasticity 

2»34i Plasticity 

Bnp ineering plasticity is concerned mainly \/ith the 
practical o. j )li cations of Lhe theory of plasticity. Y/hen a 
solid ID subjected TO external forces, it usually deforms. 

If the body does not regain its original shape after the 
external foi-c^.'S have been removed, then the body is said to 
have underjiiK jilastic def orrflation, or the plastic flov/ of the 
material in '.jie body has occured. Thus plasticity deals with 
the hehavLOur of deformable solids, v/hich undergo permanent 
deformation under the action of external influences. Only 
time LndcpCii.d oil t defonnations are considered here, 

2.3.2. Limit Analys is 

Plasiic analysis and limit analysis are synonymous 
terms. Plastic analysis is a branch of structural analysis 
dealing 'Vitli one structures loaded into the plastic range. 
Only those structures whose relationship between the generali- 
sed stresDos and generalised strains can be idealized by 
either elast lo-perfe ctly plastic or rrgid-perfe ctly plastic 
behaviour arc considered here . 
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The object of plastic limit analysis is the calculation 
of collapse loads, at which the structures continue to deform 
while the loo-ds remain constant, in other vYords unrestricted 
plastic flow occurs in the structure. In this thesis only 
the limit approach to design and analysis is adopted, i.e., 
the conditions at the point of impending collapse are consi- 
dered in determining the collapse loads. As the analysis 
and design are based on limit approach, they are termed as 
limit analysis and limit design respectively. The object of 
limit analysis 'lay also be stated as determining fclie safety 
against collepse. limit anal^/cis is based on three fundamen- 
tal theorems vis,, the lov^er bound (static) theorem, the 
upper hound (hinematic) theorem and tne uniqueness bneorem. 

In the limit analysis, the safety factor usually known as 
load factor is the ratio of collapse load to the design or 
working load. Although the calculated safety factor in limit 
anolysiG is very important, the disadvantage of limit analysis 
IS that no information is available on the behaviour of a 
structure at or below workrng loads. The ultimate object of the 
limit analysis is to find safety factor against collapse of 
a structure, but this may not be "easiblo in all the cases. 

In such cases, upper and lower bounds on the safety factor 
should be determined* In t’lis thesis both analytical and 
numerical mcthod^s are considered v^/rth reference to reinforced 
concrete cylindrical shell roofs. Emphasis is made on the 
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determinatioii of lower 'bounds to the collapse load using 
lower bound theorem. In limit analysis the structure is 
fully specified and the maximum safe load according to some 
criterion is to be determined. 

2.0.3 Terminology 

In this thesis the terminology and notation closely 
follows save and llassonnet (22), 

Perfectly Plastic Material : 

In tills moterial, strains are independent of time and 
it does not cidii'bit strain hardening properties. It is 
capable of flo\/ing indefinitely under constant stress v/hen 
the stress reaches a certain magnitude known as yield stress. 

A material characterized by those properties is termed perfe- 
ctly plastic, 

Elastic-Plastic and Eigid-P las tic Bodies : 

Both ti.e idealizations of elastic perfectly plastic 
and rigid perfectly plastic mateiial are appropriate for the 
purposes of limit analysis, fhese are shovm in Pigs. 2.2a 
and b respectively. When the elast ic-perfe ctly plastic 
idealization ic used, the limit state corresponds to the 
incipient plastic flow, .'/hen using the rigid-perfe ctly plastic 
idealization, the rigid parts of the body prevent all defor- 
mations upto the on set of unrestrained plastic flow at the 


limit load 






FIG.2.3 M-N INTERACTION 
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Collapse (limit) load : 

let a structure be subjected to a system of loads that 
are increased in proportion, starting from zero, When the 
load reaches c particular magnitude, the deformations of the 
structure increase beyond all bounds under a constant value 
of load intensity, I'he corresponding value of the load is 
defined as collapse or limit load. 

Proportional Losding : 

If all i/he forces acting on a structure can be expressed 
in proportion to one parameter, then the loading is said to 
be proportional, 

2.3,4 Physical Conditions 

Those are also known as the constitutive equations 
or stress-strain relations in the theory oC perfectly plastic 
materials. They consist of t\70 parts, the yield condition and 
flov/ rule , 

Yield Condition : 

A yield condition is a hypothesis concerning the limit 
of elastici-ty under any possible combination of generalized 
stresses. It is a function the t relates the generalized 
stresses at j^icld and is given by 

P(13) = , 


« • • 


(2.19) 



22 


where Q u the generalized stress vector, and is a material 
constant, vhiicli characterizes the yield. 

The sign of the function p adopted is such that no 
yield occurs \,iien 

< c^ . ... (2,2o) 

The yield occurs only v/hen equation (2.i9), is sabisfied and 
the comhinabions of stresses corresponding to 

nm > 0 ^ ... ( 2 . 21 ) 

are inadmiss ih le . 

The point, curve or surface corresponding to equation 
(2.19) IS knov'/n as the yield point, the yield curve or the 
yield surface respectively m one-two-and n - (n ^ 3) dimensional 
stress space represented hy the components of the stress vector 
Q. If n > 3, the surface is a hyper-surface. The yield surface 
can he shown to he alv/ays convex. 

Plow Rule ; 

The riow rule e^qiresses the ratio of plastic strain 
components during yielding. The flow rule is assumed m such 
way that the genejalizied strain vector q is given hy, 

a, ^ S 1 = 1,2, . . .n ... (2.22) 

where s is azi arbitrary factor of proportionality. Is the 
yield surface is considexed to he convex and since the plastic 
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work done is akv8,ys nonnegative 

6 > 0 . 

The Eq. (2,22) iraplies bhat the strain vector is nonnal to the 
yield surface '^iven by Eq. (2.19). Hence the flovv rule is 
sometimes referred to as the * normality flow rule’, 

2.3.5 Theorems of Limit Analysis 

Statically Admissible Stress Eieid : 

The streoses are in internal equilibrium, in equilibrium 
with thr* ai),)iiocl load P, and satisfy the static boundary 
conaitions. f ic sbresses should nov/here violate yield inequa- 
lity. A stress field of this kind ic called statically 
admissible , 

Kinematically Admissible Velocity (Displaccmont) field : 

Any cUoplucomont field compatible with the geometric 
(kinematic) houndary conditions and cerba.in continuity condi- 
tions? ly said to be a kinematically admissible displacement 
field. 

The rigid plastic model (Eig. 2.2b) and elastic-plastic 
model (Eig. 2.2a) yield the same load capacity values. 

Lov/erbound j^Statical) Theorem i 

Any load corresponding to a statically admissible 
stress field, is smaller than or at most equal to the limit 
load Pq 
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^ < ^0 • 

Upperbound (i.iae* lO Lie) Theoreii : 

Any lo^'' P corresoondmG to a kinematically admissible 
displacemeni: L'lCld is la.rger than or at least equal to the 
limit load P 

0 

Uniqueness Theorem ; 

The load capacity is unique if derived from a stati- 
cally admissible stress field for which a corresponding 
kinematically admissible displacement field exists. 

This theorem may thus be regarded as a combination of 
first tv/o theorems. 

2.3,6 Yield Criteria 

4 

In 'hie literature the Tresca and von Kises criteria 
are applied to metallic plates and shells while the square 
yield criterion is applied to reinforced concrete slab element 
in certain 63Luations. A special yield condicion (22) is 
available for a general slab element v/hich is under the action 
of the tv^o b< aiding moments and the t\/isting moment. 

The yield condition (2,19) for a cylindrical shell 
element undar consideration may be expressed as 


M • • 


(2.23) 
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Yield condition is not available in the literature for rein- 
forced concrete cylindrical shell element under the action 
of all the stress resultants. In the absence of a definite 
yield condition, the following assumptions are made to arrive 
at an approximate yield criterion. 

( 1) There is only limited interaction between the stress 
resultants i.e., there is no interaction between the 
stress resultants acting in the x and ^ directions at 
any point of the shell. 

(2) The inr.iuonce of transverse shear on the yield 
cond Ltion is neglected. 

(3) The efCect of membrane shear is decoupled from other 
stress resultants and also it is assumed that the shell 
does 'lot fail due to insufficient reinforcement provided 
to t !£c diagonal tension develooed due to 

(4) The collapse of the shell is not initiated due to the 
yielding of longitudinal steel provided to resist 

in tension zone and hence the shell fails due to longitu— 


(5) 


d inal c omp r e s o i on only . 

PQ_rlurc of the shell due to budding ii/ill not occur. 


with these foregoing assumptions, 


the EiQ_* (2.23) may 


he split into the follov/ing two independent relations as 




°2 


.. (2.24) 
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and 

= °3 ... (2.2 5) 

where Co Cr, ere meteiial constants, 

(C kJ 


If at lea..t one of the t\/o xnequali ties beconiec an 
equality, then the shell elenent is in a state of fully 
plasGic condition. 


N -yield Criterion t 
c 

As a consequence of the assumption 4 listed earlier, 
the Eq. (2.24) is valid in compression zone only and is given 
hy 

= H, ... (2.26) 

C ji 

where N = is the strength of the shell m compression 

per unit lengbii. and t are compressive strength of concrete 

and thickness of shell respe cnvely . Henceforth this yield 
condition is termed as N -yield criterion. 

M.-H, Yield Criterion : 

(fi (f, 

The condition described by Eq, (2.25) gives the failure 
of a reinforced concrete section under the interaction of an 
axial force 11, and a bending moment M. 

To obtain the M-F interaction curve (22), in a yielding 
cross section it is assumed that concrete is yielding m 
compression at a constant stress - above the zero strain 
rate level, onci below that level, the concrete has cracked and 
the steel yields in tension at a constant stress a__ * 

® j 
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denotes the cross-.se ctional area of the rein- 
forcement ner uiirt lejagfch, then (Fig. 2.3) the "bending moment 

M and the aniai force H are given by 

2 


cr t 


F = 


(l - 


and 


.. (2.27) 


]'r = A,cr _a t(l_i{»). 
st sy c ^ ^ 


a t 
c 


Denoting M = '~-'a ^-^d N = cr t, 11 and h can be rev/ritten 

u e: 0 c 


as 


I'i = 211^(1 -\|) ) 


4M 


. . . (2.8) 


!'I - 11 ( U _ 1 + ijJ ) = ( y _ 1 + i|j) 


where y = 

a 1 


IS ohe ” reduced re liifor cement ratio" . 


Elimination of the parameter ^ from E(j_. (2.28) results 


in 


y(2_y) _ 2 f- (l-y) - (f!"”') I (2.29) 

' 0 ' o "4 


lU^l 


The t\/'o yield condiciOiiG dcscr''ood by ohe Sip. (2.26) 
and (2.29) ccui be represented in three dimensional stress 
space md as shown in Fj.g, 2.4a. 

For a given thiclvneso of a shell, the interaction 

curve depends on the quantixy of transverse steel. In this 
analysis the value of is assumed such tha'b = 0.236. 
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This quantity is the amount of steel required for a balanced 
rectangular sevotion according to ultimate strength design 
under pure bcniing as per I.S. Code (25). To reduce the 
complexity oj cue problem, the yield surface is linearized 
suitably in tuc plane as shov/n inl'3g.2,4b. As is 

always compressive and for the practical range of shells, 
the condition is governed by the straight lines ef or gh 
of T'lg* 2,4b smd are given as 


IT 


= b + c 


H 

( 

N 


. .. (2,30) 


Por the value of y = 0.236, the absolute values of b and c 
are evaluatoc’ as 0.833 and 2.2 5 respectively. Hence this 
yield condition lo termed as yield criterion. 

The yield criteria described by Eqs . (2.26) and (2.3o) 


are adopted in this thesis. 



CHiieTER 3 


LL3T Al'liilitSld OP CYLIliDPJCiL SHELL HOOPS 
VlTIi PEEL LOi'lGIHJLINAL BLSES 

3.1 G-eneral 

As the solutions hased upon I'he static theorem of limit 
analysis are alvvays safe, three solutions are developed for 
simply supported reinforced concrete (RC) circular cylindrical 
shell roofs subjected to gravity loading. The load is unifor- 
mly distributed. Such a loading is considered because the 
dead \(:eight (iiicludixio self weight)of shell is a major portion 
of the total load acting on the shell. The superimposed loads 
may also be ccioidered as gravity loads depending upon the 
nature of tne lood. In such situations as they can be approxi- 
mated as giMvity loads the analysis becomes simpler. 

The general problem of limit analysis, as already 
mentioned, is fie determination of the load carrying capacity 
of a structure \/hen the material just starts flowing. As it 
is difficulc to got the exact loads (identical lower and 
upper bounds) the bounds are found. Of these, the solutions 
based on th''‘ iinpcr bound are unsafe if the correct flov/ 
mecbaiiisn is not '/isualised. 

It IS v/all hnown that the upper bound solutions yield 
only collapse loads. Since the distribution of forces is not 
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knowYiy ecoi-oMic d istribu-tion of steel is not possible as tbe 
distr iln G \ Oil ol steel is preassigred, where as in the case of 
lower bnunl soluoionc distribution of stress resultants is 
knovra througnout the shell, hence steel can be proportioned 
according to t.ie requirsuients, 

Tvfo sacHs one each in Lhe long and medium category 
of shells v/ith different geometric parameters are analysed 
herein. The same are also analysed by elastic method when 
the shells oro subjected to working loads. The distribution 
of stress resultants are shown graphically in both the cases for 
the purpose ojf comparison, 

La-tor two shells are designed by both elastic a,nd the 
proposed plastic methods. The relative economy of the materials 
used IS also discussed, 

3.2 Lower Bound Solubtons 

The 3.o\Kr bound solutions arc developed using static 
theorem. According bo this theorem, any stress fielri satisfy- 
ing the spec-i-Licd houndory cojiditions, that is in equilibrium 
with the applied loads and no vdiere violates the yield 
condition is a lower bound to the actual collapse loOvd, In 
view of this statement it is possible to construct stress 
fields which provide good lower bounds for the collapse load 
of a given structure. Thus the static method consists of 
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finding ouf o. sfatics-lly admssi'ble sfnsss field covering 
the entire doinnii. of the shell. 

jniy stnti colly admissihle stress field has to satisfy 
the equations of equilibrium (2,i5) through (2.i8). These 
four equilibrium equations contain five unlcnown stress 
resultants and Hence if one of the 

unlcnowns is preassigned then the other can be evaluated by 
solving the equilibrium equations. Y/ith this m view a 
possible distribution for at collapse is assumed, and 
the remaining, stress resultants are determined. These 
stress rosultants have to satisfy the specified boundary 
conditions, it is imperative that such solutions are not 
unique. So it is always possible to improve the load 
capacity >;. itli some other distribution. Three different 
distributions for li ^ aro assumed and lower bound solutions 
are developed ror cylindrical shell roofs. Obviously for a 
valid lower bound the stress resultcints should not violate 
the yield conditions. 

Consider the shell snown in Tig. 3. la. The geometry 
of the shell and loading are symmetric about the X oocis i,e., 
at <!> = 0 and bne loni-,itudnial edges ore free at ip= (j)^, 
where P is the angle measured from the crown to any point 
on the shell cross section and is the scmicentral angle. 
The cross section of the shell is divided into tv/o zones, 
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VIZ., the top compression zone for 0 <.<!>< p and the bottom, 
tension zone for p < ‘i’ < v/here p is an angle defining the 
location of t ix* utcil anis . In Iig. 3, la these values are 
shown. The Digress resultants have different expressions in 
the two zones hut have fco satisfy continuity conditions at 
the region boundary of the two zones. 

3.3 SyrnmetryjBoundary and Continuity Conditions 


The stress resultants have to satisfy the following 
conditions ; 


(1) 


Symmetry Conditions : 




tA 


dU dM 

=0 at d = 0 

6 d 0 d 


... (3*l) 


(ii) ConfcLnuity Conditions : 

The stress resultants IT Ih and M and the 
ol'l &1.1 

derivatives — and are continuous at <;>= 6 . 

a 9 0 9 


. .. (3.2) 


(ill) Boluidary Conditions : 


xb b 


M 


0 at 0 = Oq . 


... (3.3) 


The olielL is considered as a simply supported at 
the ends. At the simoly supported ends IJ^ = 0 at x = 0 


and X = 1 
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3.4 Lower Bound Solution 1 

In oLtaiiijuig this solution 11 distribution in blie 
transverse Cirocbion is assumed as uniform across the cross 
section both in tension as /ell as coiupression zones (Pig* 3.1b). 
As IS uniforiii across the section, ib is independent of 41 . 

This distribution can be expected in long shells if the shell 
fails as a beon of thin curved cross section. lYUen 11^- yield 
criterion governs the failure (collapse), then the collapse 
mode of a long shell is similar to that of a beam. 

The magnibude and longitudinal variation of are 

evaluated froi-i the sbatic equilibrium conditions of the shell, 

considered as r simply supported beam subjected to uniform 

gravity loading of intensity p. Then the stress resultants 

N .IT . Q li'd li, are obtained using the eo^uilibnum 
x(J)’ ’ ^4 b 

equations (2*15) through (2.18). 

3.4.1 Stress Resultants 


Evaluation cX j.T^ : 


If and represents the magnitudes of in the 

compression and tension zones respectively (Pig. 

• # « (3 #4:3.) 


11^ = - 1^}[A 


N = IT. 
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f or 0 < 4> < p » 
for p £. 'J’ 1, ‘J’o 


... (3.4b) 
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The e q_i' ilihr lam condition that algebraic sum of the 
forces at any .mction is ec^ual to zero i.e., I IT ds = 0 
and taking coii jression e,s negative and tension as positive 
yields 


2 E-( <>Q - p ) - 2 E p = 0 } 


or. 


I'l- 


^XA ^ 

103 = tt:-T) 


... (3.5) 


le t a be the lover arm between the tensile and compre- 
ssive forces diri is given by 

E( sin p - p sin 

^ ( g - p ) 


a = 


The G;:ternal bending moment Li for simply supported 
snell subjected to uniform gravity loading at a section a 
distance s is 


M = p R $_(x) (li - x) 


IS 


The moment of resistance at a section a distance x 


= 2 R p a . 


Equating the external bending moment M and moment of resistance 


M it IS found 
r 


®XA = XaTT U - m. 


• * • ^ 3 • S 3/^ 



J'rora Eqs . (3.5) and (3.6a) 


“XB = kn (1 - ^) ... (3.6T)) 

where 

a = sin p -> p sin 

Eqs. (3,4a), (3.4b), (3.6a) and (3.6b) give 


n 

^ ^ ^ (L - x) f or 0 < <> < p , 

... (3.7) 

0 cf) 

^‘x = ^ (L - x) for p < <) < 

... (3.8) 

Evaluatioi' oC I. , ; 



The N ^ js ev'xiuated in different zones as Collo'v/s : 

X^) 

The E'i . (2.15) may bo rearranged as 




. . . (3 *9 ) 


(i) Eo^O < => < p 

hif Lcrenc lading il ^ given by Eq. (3.7) with respect to 
and substituting m Eq . (3.9), one gets 
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Eq. (3.10) given on integration with respect to ‘P as 

P P 

= ""Ea* ^ ^ ^ (i* *" t- D(x) , ■ « « (3 . ll) 

when D(x) ic Pun ct ion oP x only. The value of D(x) is 
obtained by using the condition I = 0 at 4i = 0 . Thus the 

Xv 

value of P(x) is zero. 

( ll) Por p < <> < 

PifPo rent in ting given by Eq, (3.8) vrith respect 

to X and oubs'uJtuting in Eq. (3.9) yields 


^*‘x^ 

'~dV 


P <!>^ 


(1 - 2x) . 


... ( 3 » l2 ) 


On integrotior \'/ith respect to <}' Eq, (3.i2) gives 


P 


x<;> 


= - p(l - 2x)<i> + E(x), 


... (3.13) 


where E(x) is a funcbion of x only. E(x) is determined from 


the condition li,^. = 0 at b = <!>^» Thus 


E(x) = 


P 

^aT 


(L - 2x) . 


Evaluation op i!^ : 




IS Gveluatcd in the two zones as follows i 


Eq. (2,i6) IS rewritten as 


d 


I'l 4 - 1 = p, — 2p R cos 

5 5x 


,2 


« «. ( 3 f l4 ) 
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( i) ror 0 < <) < p 

(3.14:) it IS obtained that 

'!>o( <^0 “ ^ ) - 2 pR cos (<) . 

... (3.15) 

The solutiC'i lo the differential equation (3.l5) is given 
by 

nT? 

= A ^ sin cos « - pR ^ sin (j) t ) > 

... ( o • i 6 ) 

where and are constcnts of integration. On difleren-. 
tiation one go lo 

dih 

— - = A^ cos (> - sin <j-pR(sin <}>+(}) cos <f> ) . 

... (3.1?) 


Rroii Bqs . (3.7) and 


il, 


2 

3 IT^ 

riT 

rl 


P -5 

Sl, 


(ii) Ror p < 7 < 7 


Substitution of Bn. (3.8) in Bq. (3.l4) yields 


2 

an, 

d4> 


pR 

“ <!>o P - 2p R cos <}> 


... ( 3 , l8 ) 


The solution to the differential equation is 


N 


Aq sin (p + Bp cos <;> 


pR 


^ V 'r ^ ^ 


p B - pR 7 sin p , 


... (3.19) 


where Ag and Bg are constants of integration. On differen- 
tiation vhaich yields 
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dl'i^ 

= --^o cos <j) _ Bp sin <j) - pR(sin <;) + ‘l> cos <J>) , 

... (3.2o) 

The cons tents B^, Ag and Bg are determined from 
the conditions -tiveii hy Eq.s . (3.i) tlirougn (3.3) i.e., 

(i) =0 at = 0, (ii) N({, and — ^ are continuous at 

= ^ , (ill) I'T^ = 0 at <) = Thus 

A, = 0, 

pR 

= ^11 = r " ^ 0 ^^ ^ 0 " *^0 p )’ 

, 2 
Ag = ..gg pR = sin 

pR 

Br = 13gg pR = — P cos 


Evaluo.tion oi Ii^ : 


irj Gvoiaated as foliov/s : 

The Bp. (2.1.7) majr he reeirranged as 

0 2 

= —Pa R — pR COS . 

6« 

( i) Bor 0 < <!> < p 

Prom Eqs. (3.l6) and (3.2i) one gets 


... (3 .2l) 


I o o p jEu 

pR cos Cf + pR (<!> sin d - cos <>) _ — ^^( <^^ _ ^ ) . 


b4> 


a 0 ^ o 
... (3.22) 



41 


On integration onc.e it is found 


Oivi g 2 

— ~ = -^ 1 -, pR Sin pR <{> cos <> 

6 (j) 11 


Second integration yxQldi 


)0 + 0 . 


... (3 .S3) 


2 2 'T ^ “ 

pR coo <) - pE ^ )+cos + <f> sin <i> +C^ ^ +0, 




... (3.24) 


where and Op aire cons bunts of integration, 
(ii) Por p < ^' < ij) 


Sahc bite t ion of Eq. (3.19) in Eq. (3»2i) yields 


S^'^A 2 ? PB-^ 

— = -Agg pE sin (t) ~ Bgg pE'' cos p + •— p 

5 0 

2 

+ pE" ( 0 sin 0 - cos 5 ) . 


... (3.2 5) 


first and second integrations of Eq. (3.2 5) give 


2 2 pE'^ 

= Agg pE cos 0 - Bgg pE sino + ~ 5^ ^ 5 


pE 0 cos 0 + C- , 


... (3.26) 


2 2 pR 2 

= Agg pR sin 0+ Bgg pR cos 0 + 5^ p $ 


pR ( cos 4) + <j) sxn + j (3 #27) 


respectively. and are constants of integration 
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The constants C^j G 2 * ^3 ^4 obtained using 

the conditions given by l4qs . (3*l) through (3«3) i,e*, 

511 6 M 

/ 1 \ _-«ll _ 0 at <l» t= 0 , ( 11 ) 1 '^* s-iiLd are continuous at <}> = Pj 

(ill) = 0 at 4 > = < 1 >Q. Thus 


C, = 0 , 


C. 


c, = 


^22 


pR = 


pR 

3. 




-p) (1+ “) + a cos ^ 


2 2 
■^33 “ ” "sT" "^0 ^ * 


C 4 = pR = 


2 pR' 


a 


a 


cos p ( 1 - 


Evaluation oT 0 , : 


Q(jj ID evaluated as Pollows * 
Eq. (2.18) IS rewritten as 

1 t 

V.TT * 


IromEqs. (3.23), (3.26) ana (3.28), at as lound 


Q = -B_ pR sin (f) - pR 
4* 1 -L 


a 


( 4) ) (j) + 4' cos 4 


... (3.28) 


for 0 < <)> 4. 

... (3.29) 


and 


Q ^ = A22 pR 4 > 


^ ®22 4 + pR (-“-^ - ^ cos 4>)--Cg3 pR 


for p4<}><4>Q« (3 *30) 
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Tlie s'bress resultants after substitution of the 
constants of rntegration boil dov/n to : 


(i) lor 0 < ^ < 


^ *^0 / A 


(x) (1 ~ X) 




p <J> 


x<{> “ ^a ^ 'o 


° iK - Cl - 2x) 


... ( 3 . 31 ) 


... ( 3 . 32 ) 


= 


pR 

a 


a sin 4 - eos <{i^ 


cos ^ ) cos ({> 




... ( 3 . 33 ) 


Ivl. = ~ 


££ 

' a 


2 


(cos <b -v 4^ sin (> - cos 4,^) - <))q(P cos 


^ cos p ) cos 4 > + 2“ 




Q , 


pR 

a 


4) cos 4> + 4 > cos 


+ ‘l>o - 


... ( 3 . 34 ) 


I in <(> 


... ( 3 , 35 ) 


(11) lor ^ < « 5 . 

PC'-) 


P '>0 


^x4) " ^ 2a " '0 


^ ( 4> - 4> ) (1 ” 


pR 


= - a 


a <!> sin 4> <;)q sin ^ sin 


... ( 3 . 36 ) 
... ( 3 . 37 ) 

^ + 4,^ p(l-cos cos <}> ) 
... ( 3 . 38 ) 
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= 


*a~ 


'(cos cj) + <^ sin 4' — cos 4>) _ <!> sms sin 4 

0*^ o ^ 


4’q P cos cos 4> - - 2 ) 


0 • m ( 3 • 39 ) 


e£ 

a 


2 

0 4) cos 4) Gin p cos ^ 4^^ p cos sin 4> 


+ ■fo |3 (■>„ - ■» 


• • ♦ (3 *40 ) 


The e:cpressions for the stress resultants , LI 
and Qji, indicoitc that they do not vary with x. 

3.4*2 Detenainat ion of limit load 


The liuiit load is computed from the tv/o independent 
yield conditions and the minimum of the two yalues is taken 
as the collapse load. Let p^^ and represent the t\7o 
lower hounds to the collapse load found from and 
yield criteria respectively. They are evaluated as follows : 


(i) Due to IT -yield criterion : Prom Eq. (3.3l) it can 
c 

easily he found that N is maximum at the centre of span. 
Thus 


H 


X .'IJ; 


P 

8aE, 




) 1.“ 


. .. (3.4i) 


Yield condition (2.26) g,ives 

8aK IT 

c 

p _ 

<i>o (<^n - P) ^ 


• ft 


(3.42) 
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Then p^, the limit load m nondimensional form is obtained 

after dividing by a as 

c 






(3.43) 


(ii) Due to yield criterion ; The failure of a 

cylindrical shell occurs in the transverse direction where 

¥ 

IS maximum. In a simply supported cylindrical shell v/ith free 
lonf^itudinol ed^es, is maximum at its crown. Botn 
and do noc vary in si^n ior the entire shell and li ^ is 
alv/ays compiessive . Since and are independent of x, it 
IS understood that if once yielding is initiated at a point 
then the enlire shell yields along the crown. Let 17 and 

Y ^ 

represent the values of M<(, and at crovm. Then 
Eqs. (3.33) nnd (3.34) yield 

" r L *^0 °°° *^0 "■ CO® P) + <^0^ <>0 ” ^ ^ J ’ 

... (3.44) 

P 

M. = - - a(i cos p^) - cos cos p) 

- 2° ( <!>o " P ^ ‘J’o ^ * ... (3 .45) 


Let pp denotes the nondimensional load, pg is found 
from the yield condition given by Eq, (2.30) using the 
Eqs. (3.44) and (3.4=5). Thus 
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Po ~ 


"12 ba 




* • • ^3 * 46 ) 


where 


\ = 4(f) “ 


a ( 1-cos _ (t^(p cos _ <j,^ cos p ) 


- ^ (i-0 - + 2) 


+ o(|) 


* cos i 


♦ o oos p - ^ ) 


Failure of the shell occurs either in the longitudinal 
direction due to II -yield criterion or in the transverse 

o 

direction due to yield criterion whichever is critical. 

Denoting the lo\7er bound in nondimensional form Pq 
IS equal to the ruinimum of the two values p. and p^. 

3.5 Lower Bound Solution 2 

In che solution i, II is tal-cen as uniform across the 
cross section, \tiich requires uniform distribution for steel 
in longitrdmol direction in the tension zone at collapse. 

But the elastic analysis reveals that distribution is 
almost linear in the tension zone at v/orkuig loads. Therefore, 
the assumed uiiiform distribution for IJ ^ may not be satisfa- 
ctory. In view of this the distribution for IT is modified 
in the solution 2 as shown in the Fig. 3,ic. In tcis, the 
B IS uniform across the cross section in the compression 

IHZ 

zone and it varies linearly with depth z' from the neutral 
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axis in the tension zone. The stress resultants are 
evaluated on similar lines, 

3,5.1 Stress Resultants j 

Evalua.tion of IT, : 

let he the vertical height of the centre of total 
compression from the neutral axis (Rig. 3.1c). It is same 
as the centroid of the circular arc of the shell above the 
neutral axis . Thus 

- .sin S , 

h^ = « cos p ) . 

Let h^ and hg he the depths of compression and tension 
zones (Pig. 3, la) respectively which are given hy 

h ^ = R( 1 - cos p ) , 
ho = R( cos a _ cos P ) . 

If z’ IS the vertical distance from the neutral axis 

to any point across the section which is given as 

a* = R( 1 - cos <!)) for 0 < <J> < P > 

z' = R( cos <J> - cos p ) for p <_ (> < . 

The maximum value of m the tension zone and the 

uniform value of IT in the compression zone are denoted hy 

yL 

R and respectively (Pig. 3.1c). Then the variation of 
P 

R in the tv/o zones is expressed as 
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1 


X 


- I'l. 


,CA 


f oi’ 0 < < p 5 . . . ( 3 ,4 7a) 


N = 

X 


fNj r- I 

j?. :i 

hg 


cos ^ _ cos 7 ) 
( cos p - cos (J ) 


fox p < ... (3.47b) 


The St uic onuilabriufli condition i.e., £ 11 ds = 0 


at any section o. yields 
2 


IT E d<J) + 2 / ° IT E d<!’ = 0 

311 _ IX 


... (3.48) 


Subst itut iii{i from Sqs . (3.4 7a) and (3,47b) and integrating, 
the following ns obtained 


C( ^ )cos p - sin <j) + sin p 

2E1'T.v-a & “+ — = 0, 

^ ( cos p - cos <{> ) 


‘XI 

which yields 


N = 
P 


P ( cos p - cos i>^) 

[3 ( <(>Q - ^ ) cos p - sin 4-0 + sin p 3 


. ... (3.49) 


Tahing moments about the neutral axis, the moment of 


resistance 11^ at a section a distance x is given by 


k 




d E d<?> z’+ 2 j IT E db s’ . 


... (3 , 5D) 


Substituting from Bqs . (3,47a) and (3.47b) and performing 

ui. 

integration, Ii is evaluated as 


Jg = 2 ii^p E a 


’ 1 ’ 



where 


1 , sui a, sm 2* sin ap 

h = a L~2'-- (1 -T^) + -- X-- + —4— 


(sin p sin <1) ) - cos p sin 'J' 


= ( - p ) cos p - sin <) + sin p 


The onternal hending, moment li at a section a distance 


X, under uniToJ’ii gravity loading p is 


;E <>q(x) (I - x) . 


Equati 4 i,_ the external hcnclmg moment M and moment of 
resistance IT one gets 


p c>q(x) (L - x) 
"2 a^ p rr 


(3.51) 


■Using hqs. (3 , (3 .49) and (3.5l), li,. in the tension zone 


is found. Thus 




(cos p - cos (;> ) (x) (L - x) . 


... (3.52) 


SuhstitutiOi- ofEqs. (3.5l) and (3.52) in (3.47a) and (3.47b) 


respectively yield 

P <>Q(n) (L - x) 

H/ a*k ■TM i i.t* ' newfenri. » 

^ 2a p R 


for 0 < <!> < p j •• • (3.53) 


p 

Yl _ ^oos p cos <))(x)(L - x) for p < <> < 

^ 2a^aR 
1 


... (3.54) 
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Evaluat lOii 


0± iJ 




li , 10 evaluated as Xollows i 
x<p 


Ec! . (-. 15) gives 


3 • 
6 <i> 


6N 


= - R 


X 




... ( 3 . 55 ) 


(i) Ror 0 < ^ ^ 

Eqs . (''. 53) and (3.55) yield 




... (3.56) 


On integration \;itn resped to 4) it is found 


x-!) 


l-e~r 


... (3.57) 


wliere D(o:) lo a function of x only. D(x) is J ound by using 
the condition 0 at <>=0. Thus the value of D(x) is 

zero . 


( ii) For p <; (j) < 

Elis, (3.51) and (3.55) give 
fill p (i> 

(cos 6 - cos 3 ) (L - 2 x)... (3.58) 
e 3 Eaa^ ^ ^ 

On integration \/ith respect to 3 one gets 

-{3 cos p r- sin (|))(l-2x)+E(x) , ... (5.59) 
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where S(x) is e- function of x only. S(x) is determined from 
the condition . = 0 at <j> = <> . Thus 

X (p ^ 


E(x) = “py™ (Oq cos P - sin <J)^)(L-2x) 
1 


... (3.60) 


Evaluation of : 


IS ova-luated as follov/s J 
Eq. (2.16) IS rewritten as 

2 T - 2 __ 

e Ha 2 ^ ^ 

I'f^ -j. ^ 2pll cos (J) 

6 x 


... ( 3 * 6 1) 


( i) Eor 0 < 6 < p 


Substitution of Eq. (3.53) in Eq. (3.6l) yields 


6 Ha pE 

* r/ = 


~ 2pR cos <) « 


... ( 3 . 62 ) 


The differential equation has a solution of the form 

= A, Gin<^+ B. cos <{. + - pE ^sin c? , ... (3.63) 

<>1 1 P 

Where A end B are consbants of integration. On differentia- 

1 

txon one gets 


^ ~ A ^os « - B. sin <> - pE(sin <- + 4 > cos ^) . 

60 ^ 1 1 ... (3.64) 


(ii) For p < 0 < 0, 


-Otis. ( 3 , 5 t) a»a 43,611 give 






5 8.17 
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d i'u 

H + - 

0 tj) 


pR d 


aa 


I 

~ (cos p - cos <)) - 2pR coS({) . 


... (3.65) 


The solutiov- oC the differential eq_uation is 


H + Ag sin ({) + Bg ^ “ P®- 


+ (J) sill <J) 


cos 


3,0. 


2aa. 


<!> sin <> 


... (3 .66) 


■where Ag a,nd Bg are constants of integration. On differen- 
tiation Eq. (3.66) gives 


dN 


d 0 


^ = A^ cos 0 — Bp sin 0 - pR(sin (p + <p cos ^(l - 2a‘a“^ * 


... (3.67) 

The constants A,, B., Ap and Bp are determined from 

1 1 ^ 5P 

the condibions (3.i) through (3.3) i.e., (i) = 0 at 


dll 


^ = Of (if) llj, continuous at (^= p, (m) = 0 

at 0= The constants are rewritten as A^ = A^^ pR, 

B = S- 1-1 P ®-7 -^g " -^22 P^ ®2 “ ®22 P^* where 
1 11 


Ai^ = 0 , 


B =z (p tern 6 


4 cos B ^ A 


- nos ^ - 3 


• ( sin ^ 6 q + # ) » 
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^22 - (sxn 8^ _ 2^ ) + ^° , 


Egg = tan 


1 - + air- (t-,-F^ 2P - 2 2? ) 


a^P 4aa^ ^sin 


Evaluation ox 1,+, : 


li(jj i& determined as given below 
Eq . (2.i5') may be rewritten as 




E. R ^ pR cos <) 


... (3.68) 


( i) For 0 < ^ < P 

•Kixiw- 


Suustitution of Eq. (3.63) in Eq. (3.68) yields 


1. 


6 ri 


— dR"'’ cos 6 - pR (’g~^ *’ *!’ sin ^ )-pR coo 6 


66 


11 


■# 


. . . (3 .69) 


On first integration it is found that 


Si'ia, p 2 '^n'^ 

pr sin 6 - pR (*~°-j + 6 cos 6 ) + C^. 


.. (3.70) 


Second integro-tion gives 


pR^ cos 6~ '*' P°® 6 + 6 sin 6 )+C^ 6 + Cg, 

... (3.71) 


where 0. and Gp are constants oE integration. 
1 ^ 
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Cl, = 0, 

<^n 


<>r 


C33 = - aa- ’ 

C44 = cos <^Q + a, ” aa. 


cos 4 >q + (1 + 2 ^ 


+ (sin ^ - P cos 
•jio / . * 2^ _„ . *^o / - . 2. 


(d t 0 ^^) coS^-K (24^ + 0 ^ 3 ^ + C 44 . 


■'22 2aai ' ■ 


Evaluation of Q^j, : 


is determined as follows i 


<l> 

Eq, (2,18) IS rev^rrltten as 
1 

'rW * 


... (3.75) 


EromEqs. ( 3 , 70 ) , ( 3 .73) and (3.75), it is found 


for 0 < ‘I’ d 


. . , (3 .76) 


and 


Q ^ d,Q Q P ^ cos 


,j) » Egg (j) + pR 


r "COS p 


aa 


<t> 


1 


2aa 


1 


(sm 4’— 41 cos 4^ ) - 4 oos 4> 

for p d 4) d. 4o 


4 - C 33 pR 
... (3.77) 
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The stress resultonts are sunmarrsed below : 


( i) For 0 < (^ < p 




(x) (L _ x) 


... (3.78) 


xffi Tp 


i|> (1 - 2x) 


... (3.79) 


~ pF j cos ^ ^sin <j) 


... ( 3 . 80 ) 


_o r ^0 <!> 

cos <t> - cos (()+•<() sin <j) + - Gr 


... (3.8 ]_) 


= -.pE sin <)>+<}> cos <|) -f- ^ 


... (3. 82) 


(ii) For p < <)> < 4» 


*x “ T'faa" R 


( cos p - cos )(x)(l 


... (3.83) 


Taa’^ L p + sin 


sin (L~ 2 x) 


... (3.84) 


iin <|) t£~" •(cos p 


<i> sin 4> 


) _ A p sin 4 » 


- loo cos ({) 


... (3.85) 


It = -pi" cos (j) + 4 . sin 


(2 cos <{> + 4 ) sin <{) 


+ 4 *" cos p) - Agg sin 4> Bgg cos 41 - 


... (3.86) 
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Q , = -pR 


Y 

4> cos 4) _ 2 ^^ ( 4) cos 4 ) - sin 4 >+ 2<J> cos p )- 


^22 4> - C 


22 


'33 


(3.87) 


3. 5.2 Re beri'iination of Limit Load 


The limit laad is determined by using the yield 
conditions. Let and represent the two lower bounds 
found from and yield criteria respectively. They 

are evaluated as follows ; 

( i) Due to li — yield criterion ; The maximum value of II _ in 

O Ihu 

T 

compression is determined from Eq_. (3.78) at x = ■^ . Thus 

A 

P 

N .. = ^ . ... (3.88) 


X max 8a * 


The yield condition (2,26) gives 


So pR H 

Pli - "L'P'" 

^0 


... (3 .89 ) 


Denoting the nondimensional load as p^, I'c is gi\re 


3 ^ as 


P 


1 " «y|)hf) 


(3.90) 


(ii) Due to il^ yield criterion : Failure of the shell 
occurs ob the conra in the transverse direction (see section 
3.4.2). Let I, and represent the values of and 
at crown i.e.,at 4 = 0. Eq.s. (3.80) and (3.8l) yield 
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— — pR (l — “ ^22^* ••• (3»92) 

Let pg represents the load xn nondimensronal form. Substi- 
tuting and in the yield condition given by Eq_. (2,30) 
it is found 

_ £i2 _ , . . ^ 

■^2 f “R 2 P ‘j’n n * 

j_ 4(|) ll-Ogg^^P 4- o(f)(B^^ + 

... (3.93) 

The minimum of the tv^^o values p^ and Pg lo denoted by 
p^, which is the actual lower bound. 

3.6 Lower Bound Solution 3 

In the lower bound solutions i and 2, is Lohen as 

uniform in the compression zone. This uniCorm distribution 

for R IS oOGOible only if the shell fails as a be.ui of 

curved cross section. Depending on the geometric parameters 

of the shell, ib is possible the failure of the shell in the 

transverse direction due to yield criterion nay precede 

the I -yield criterion. If this failure occurs, then the 
0 

distribution of may not be uniform across the section in 
the compression zone. Logically the distribution of in 
the compre^ion zone is assumed as parabolic, varying with 
vertical distance from the neutral axis. In tnc parabolic 
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variation, fclie maximum value of N in the transverse direction 
occurs at an angle (Pig. 3,ia). This angle p ^ is kept 
as variable and its value is determined for a shell of given 
geometry sucxi tiiat it gives optimal value for the ultimate 
load. The 0 is oribution for II in the tension zone is taken 
as linear venation as in the lov/er bound solution 2. This 
IS shown in Pio. 3. id. 

3,6.1 Stress Besultants 

Evaluation oT * 

I'l jc determined as described below : 

1 end il represents the maximum values of IT in the 
IP" X 

compression and tension zones respectively as shown in the 
Pig. 3. Id. As described earlier the variation of in the 
transverse direction is given by 



for 0^4* ••• (3 .94 ) 


for p < 4> .«♦ (3.9 5) 


where z' is the vertical distance from neutral axis defining 
1 

the position of the maximum compression and hg is depth of 
the tension zone^ Thus 


z’ =. R (cos cos ^ ) for o < <}) £ ^ 5 
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z’ = E{ cos ^ - cos for ^ < <}• £. ♦qJ 

z^ = E( cos p ^ - cos p ) , 

Eg = E{ cos p - cos ♦q) • 

These are m^rcated in the Pigs. 3.1a and 3. Id. On substi- 
tution of tue a;oove in Eqs. (3.94) and (3.95) it is obtained 


IT 


X 


—S' g (2 cas p cos p 

( cbs p ^ aos p ) 

- cos p ) for 0 < <!> < P 


cos 0) ( cos <l> 
... (3.96) 


and 

L'l ( cos 6 - cos (j)) 

II - 

^ ( cos p - cOS i^q) 


forp £,<)> S'i’o* 


... (3.97) 


The atobic equilibriun condition i.e.,i: ds ^ 0 

at any se colon of tlie shell yields 


2 / y i: d <l> + 0 / 1; Ed 4' =- 0 . 


4>. 


0 


o 


j.^ 4 'rro Qf^') and (3.97) and on carrym^^ out the 
Substituting ncis^ (3*9b; ana ; 

integration it is found 


-2 r ^ ^ J ^ ^ cos p) + g cos 2p 

( cos P ^-cos p ) ^ 

JT E 


sin 


_„-J2 L(‘*’o"^^ ^ 

( cos p “ COS (J)q) 


- sin + sin p " 0, 



6l 


which gives 

iT^ ( cos p - cos 
^ a^( cos ^ ^ - cos p )2 


• * » (3 *98 ) 


w'here 


a = 
1 


(‘J>Q “ ^ ) COS p - Sin (fi Q + sin p 

/ V ^ 

con ^ sin p - p cos p ) + ^ cos 2p 


sin 2p 

"■“T“ 


Tne moment or resistance at any section is found by taking 
moments aboiil: the neutral axis. Thus 

P 

Vi = 2 f IT R d <!> z + 2 i IT R d <}) z’ . 

^ 0 P 


Substituting Ilq.s. (3.96) and (S.gV) and performing integration 
one gets 


11 


2iT^E ag 


( cos p ~ cos p ) 


’2 ’ 


... (3.99) 


v/here 

ag = p cos p ^ (2+cos 2p ) + sin 2p cos p - f cos p 

- I sin p - !■ cos p cos 2p + C ( “IJq - P ) (^+2 cos 2p ) 
+ sin 2<(>q - 8 cos p sin + 3 sin 2p[] . 


The exterxiil bending moment M at any section x, subjected 
to uniform g cavity loading is 
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I'l = pR ♦q(x) (L - x) . 

Equating Bqs ♦ (3.99) and (3.100) one gets 

P K 2 

^^q ^ ^T1 (cos p ^ - cos ^) (x) (L - x) 

2 

Using Eqs. (3.96), (3.9V), (3.98) and (3.i0l) 
expressed as 

P 4>q 

(3 cos ^ ^ - cos ^ - cos ,j) (cos <l> 


for 0 < ^ 


and 


H 


X 


P ({i ( c jS p - COG ({) ) 

4x)4L-x) f or p < 4> 

1 


Evaluation of i 


11 ^ iG determined es described "belovi/ ; 

X 9 

Eq. (3.15) IS rev/ribtenas 


&4> 


R 


di 

6x 


X 


( i) Eor 0 < <<> < p 


Substitution of Sq, (3,102) in Eq. (3 


5 f ~ 


cos 


^ - cos^- 


cos 4) ) ( cos 


(3.100) 

. • » * ( 3 . 10 1) 

, U can be 

* X 

- cos p ) (x) (L-x) 
... ( 3 . 102 ) 

~ ‘i>o**** (3. 10 3) 


(3.104) 

104) yields 
(j) - cos p ) (L-2x) . 
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On inte gr otxn^ ^7i-th respect to (j, 


N 


P r 


xi}) 2ar 


o 


cos ^ ^(sin 4) -^ <}> cos p) + ^ cos 2p 

(1 - 2x) + D(x) , ... (3.105) 


Gin 2(f> 

-* —r'™- 


where D(x) is a function of x only. 

( li) For p 5. ‘I’o 

Bqs. (3,103) end (3.104) give 


3B . p <J) 

- Q- - 

6 'f *^^1^2 


(cos <p - cos 4>) ( f - 2x) 5 


v/hich on integrating v/ith respect bo <t> yields 


P 1 

xi;j ■ " 

where B(x) is 


( cos p - sm ({) ) (L - 2x) + B(x) , 
fuiLction of X only. 


(3.106) 


Evaluation of I'(j, : 

IS evaluated as follov^s : 

Eqr (2.i6) results 
S^h. P 6 B 

N* ^ *.-o- = R pR cos . ... (3.107) 

^ dx 

( i) For 0 < <(> < § 

«r w» J-®- JWf.^ me *S ^ *P»M 

* 

Frora the Eqsj (3.102) and (3.10?) ib is obtained 



64 


6 pr ,: 

If + — — (2 cos 6. - cos a ^ cos (cos <J) - cos S) 
^ d<P ^2 1 


- 2pR cos (j) 


... (3.108) 


^he diflerential equation (3,io8) lias the solution in the 
form 




dR 

= A sin <f> + B cos (f) + ™ 6 ( <f> sin 6 cos g. 

JL X ^ ~ 1 

2 cos p cos ^ *■ _ pj^ <{( sin if , 

... (3.109) 


v/hcrc A, ij'o are constjnbs of integration. On differentia- 
11 

tion the follo\;ing is found. 


91^ A dR 

= A rct, if B sin <f + 
d<f 1 1 ^a^o 


cos p ^ (sin (f + <f cos if) 


;in 2<f 

r/ 

o 


- pR (sin (f 4- (f cos <f ) . 


... (3.110) 


( ii) Ror p .i. 4> < <fc 


SuhsLitrtion of Eq, (3.103) in Eq. (3.10*7) Alves 


d^iu 


pR <f^ 
^" 1^2 


( cos p - cos (f ) - 2pR cos (f . 


... ( 0 » 111 ) 


The solution to the differential equation is 


pR pR 

= ^2 ^ + Eg cos 4 > - ~ cos p + 2 --^- if^if sin (f 


Lq J-iJ. y -f- ^ 

- pR Sin <? , 


... ( 3 . Il2 ) 
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where Ag and Eg are constants of integration. On differentia- 
tion one gets 


An cos 0 - Bn sin (J) + 


an o 


(!> ( sill <J) + <}> cos ^ ) 


lii Gin (fc + <> cos ) . 


« . * ( 3 » ll3) 


The constoncs are modified as A. = A^ . pE, = B_ oE, 

Ag “ "^2 B ~ ^2B P^* 

Evaluation of : 

Fi^ IS determined as follov/s : 

Bq. (f.lV) yields 

n-- 

S i-. 2 

-■'g* ^ -11(4, R _ pR cos 0 • ••• (3»ll4) 

3 

( i) For 0 < ^ 

From Bqs. (3.109) and (3.ll4) iL is obtained 

o 2 0 R^ 

dR'" sin 3 - B pR cos <> - (j (cos p . ‘J’ sin 3 
11 '^ 11 ^ 1 

cos 2p cos 23, ,2, ^ ^ 

COS ^ ^ cos p j 4 p};{ (p sxn 9 -cos 9 ;♦ 

... ( 3 • ll5) 

First and second integrations with respect 3 give 
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■^11 <*> ~ pR sin 4 i - -Y~ <f>Q 


6 <l> 


cos p ^(sin <> 


- () cos (^) - 2<t> cos p ^ cos P + 2 

O 

- pR^"*!) cos <J> + C. 


cos 2p sin^ 2<J‘ 


• •* ( 3 . 116 ) 


and 

2 2 

pR' sin 0 4 - pR cos <{) ~ — — (jj^ 


pR^ 


cos 2p ^2 


1 


2 cos B'J'"! 

- COG P^(2 cos 4> + (> sin <>) - (cos cos p)^> - 

• • • ( 3 • ll 7 J 


,2 


- pR ( COG 4) + (J) sin ij) + C <J> + Or 


where C .nd Cp are constants of integration, 


( 11 ) Bor ^ < <t>. 


iSuDG titution of Bq. (3.ll2) in Eq. (3.114) yields 


1 


d du 


/2 


-.2 


— = -Apo pR sill C) ~ Bgp pR cOb ^ + pR ( 


2 ^ 


d <> 


__ JD^ 


6 sin (J + 6 sill <) - oos t? ) 


°'1^2 


• • « ( 3 » 118 ) 


.c^i- 2 

On first integration it is found 


2 2 2 
r-t = Boo Oil cos (J - Bpp pR" sin (J) + pR 
09 


’22 


cos p 

4> 


a, 


12 


^ - (sin <}> - ij* QOS c;;) ) - <) cos 




+ Cg ... ( 3 » 1 19 ) 
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where le a constant oT integration. By second integration 
one gets 

= Agg pR^ Bin 

+ <{> sin 4> + 

. * • ( ^ • l2 0 ) 


6 + B^p pR cos 6 + pR 




<) COS p) - cos 6 - <> sin 6 


+ + C^, 


where 0^ is a consbcint of integration. The constants are 

2 2 2 
modified as = ^11 ^ ^2 ~ ^^2 * ^3 ~ ^33 and 



Evalua.tion of : 


IS determined as descrihed helov< : 


Eq. ( 2 , 18 ) yields 
, 1 

'^6 tt hV * 


(3.121) 


From the Bq.s, (3,ii6), (3.119) and (3.l2i) it is found 


pR 

Q = A.. pE cos - B pR sin ^ — <>q 

^ ± X ^ 


cos p ^(sin (j) 


(f- cos 7 ) - 2 $ cos cos p ^ 


cos 2p sin 2<i> 


(j) + 


12 


_ pR <}> cos (f> + pR Bor 0 < <> < p ... (3 .122) 


and 



S3 


= -^^22 ^ ~ ®22 pS sin <{) 4 . pR 


<*>Q COS p 
®' 1®'2 




Xsin $ ^ $ cos <{))«, ^ cos 4 > 


2a .a^ 
1 '-^ 


+ C 33 pll 


f or 1 1 ‘i^o* ... (3.123) 

Tho valuee "or '^11’ "^11^ “^22’ ^22* ^11^ ^82^ ^33’ 

^44 are o^o-taijaed applying^ tile conditions (o.l) through (3,3). 
They are given by 


P(x) = 0, 

= h:k p h> - ^■'-'^ ’ 
1 ^ 


Aii = 0 


B 


11 


<}> 

^ 0 

(p taxi 0 ^ + 

0 0 VC 3. 


tan p - fp cot p - 2(l+p! cot p) 




cos p (1 + ^ cot ^ ) 


2 


cos S - f . tan ‘P 
0^1 0 


cot § 


A ^ f« 

^22 = ^2 ^^2 


B 


22 


<)_ ton <>. -K Tr--““ tan p„ - — f, tan p^ , 


^0 4a ^ag -3 


0 ag 1 


'11 


0 


Cgg = Oc 


^o 


ip- 

■6- "■ 4af 


n (1 ~ ^ 0 ^ ^4 + 


9 
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n r f 

^33 = -^5 “ 4a2 ^4 » 

A a2 

<? «1>Q 

°44 = '^o - ET°5:J ^8* 

where the constants i = 1,8 are given hy 

f = 6p cos p — 3 cos sin2p - 2 sinp(l- cos 2^ ) , 

1 1 -L 

fg = Sin 2p - 2p , 

= 4 cos p cosec <)q - sin 2p - 2 ((>q - p), 

= 2p cos 2p + 8 cos p^(sin p - p cos p ) - sni 2p , 

f g = sin p - p cos p , 

fg = 2(<'Q-p} sinp +2 cos 4 |q + ( - P ) cos p - 2 co; 

o 2 

f ^ = 4(4 + p cos p ^ cos p - (2 .5 + p ) cos p , 

fg = 2 cos Oq (2 + 4 “ ^ *^0 cos p + 2 0^ sui p . 

The stress resultants ar® suiuuiarised as follows 


I'or 0 < '^ < P 


P 


I = 




28 


V 


\ (2 cos p - cos p - cos o) (cos « - cos p)(x)(L x) 

>.rt 1- . .. 


(3.124) 




= P 


Clf 


sin 2' 


cos p^(sin (^ - 0 cos p) + I cos 2p - — g j 


3 (L-2x) 


• -» « 


(3.125) 
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L V 

cos 0 - <{> sin ^ + |p(4>cos sin (j - 


- 2 cos p cos p . + 


cos 2p cos 2*^ 


) 


... (3.126) 


M. 


= pR ^^11 ^22 ““ 4> - <}> sin (j ) + 


^22 


+ pR 


2 *^0 

a. 


r 


cos p (2 cos <> + 0 sin <> ) + 


,2, cos 2a ^ cos 2^“ 

+ <> ( cos p ^ cos p « 


... (3.127) 


= -pR 


a.-, 


cos p ( sm <> - <!> cos <>) - 2<^ cos p cos p 


co.^ 2p sin 24> 

■^2 ■'* * “‘ I?'"" 


pP-CBii sin 0 + <) cos 7 - 

... (3.128) 


For 0 < <> < ’i'. 


= P 2-aVa-,l ( ^ 

1 


COB o) (3c) (L - 2x) 


... (3.129) 


w ° , 

^^x<) " -P "^^2 


sin (J - sin <1)^ + ( 7 q - (?- ) cos p 


(L - 2x) 

... (3.130) 


B. = pR 
C> ^ 


A„o sin 0 -f cos - (>sin (J + 


+ ^ S3J1 -f 2 COS p 

^^'1^2 


... ( 3 . l3 1) 
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i 


sin <(> + Bgg cos 4 + Gjg 4 + 

2 cos 4+4 sin 4 + 4^ cos p ) I 


cos 4-4 sin 4 

• • ♦ 1^ 3 • x3 2 ) 


Q* = pE 


Ao., c-ori ■* 


B 


gg sm 4 + ^33 cos 4 


4 cos p (sin 4-4 cos 4 ) 


^12 




.. (3,133) 


3.6.2 Deternination ol Limit Load 


Li.iit loo,d IS found out from tiie yield conditions. 
Let and "be tlie lower Bounds obtained from the and 
the yield criteria respectively. They are evaluated as 

follov;s : 


( 1 ) Due to ii -.yield criterion i The maximum value of F 
will occur ^t toe centre of the span of me suell in the 
lon^situdiiiej direction and at 4 = p ^ in l^he transverse 
direction. Substituting 2 : = ^ and 4 = m Lq . (3.124) it 
IS found t.io L 


•P 4 


u 


2 ^ 


X max SagE r ^ 


, ( cos ^ - cos L 


• , • ( 3 , l34 ) 


The vield conditi'ii (2.26) gives 


P-i 


Sag £ F ^ 


*^0 (cos ^ ^ - cos ^ )' 


, , » ( o , l3 5 ) 


The nondiuonsional load is evaluated as 
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P. 


Ll 


8a^ 


^(|)" (|) 


... (3.136) 


- COE 


(li) Du© uo yield criterion : Failure of the shell 

occurs a,t the G^-OTm in the transverse direction (see section 
3.4,2). let j' ^ ^ anrlM^^ represent the values of and 
at crovm i.e.,at <!> = o . The Eqs. (3.i26) and (3,127) yield 




‘<’0.1^ cof_f ^ 
^2 ^ ^ 


2 cos 



j 


... (3.137) 


I'l 


<l> c 




(3.138) 


Pg denote.© tuo loci,d :n nondime ns lonal form. Substitution of 

h Slid li in the yield condition given by Eq. (2,50) yield 
(pc pc 


Po 



b 



... (3.139) 


where 


1. 


o 


= (1 - Cgg - 1^^ - 34g_; 


2 r-" cos 6 .) 


+ c(:|:) 


T. ^ io a 4. 


a. 


_ 2 


cos p cos ^ ^ 


)]• 


Then is Lie minimum or the t./o values p^ and pg 
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3. 6,3 Optxmel Value of ^ ^ 

In uiie IDqs • (3,i36) end (3.139), p ^ is a variable 
while othej. x'^oianeters are constant. The value of is 
determined in rjuch a way that and pg a-i© optimal. The 
values of p are determined numerically, lor different 
values of ’t> values of p ^ are found for a given value of p . 
There values are tabulated in Table 3.1. 

Tabic 3,1 ~ Values of p ^ in solution 3, 

<^0 ^ 

30° I 0.1851 

0 .1091 
0.2327 
0 .1309 
0.334 5 
0.1527 
0.3818 
0 . l636 






2/3 0 , 


< CN O ^ 
l| H 11 ll 


■©• cx CQ. nr. 


i 1 

X? 

D 

O 

- 

s. 

o 

n 

u. 


o 


C9 



$ 



a ,i 

o 

O 

O 

Si 

_! 

-J 

-J 


1 

I 

1 

OO 



CM 

rr 





0001 X cl 


o ^ 
n 

O') O 


•©- crx cx o:: 


nnni x c 


oo 

d 



fc> fc> 

^ ^ 

o o. o 


f I I 

CM 00 
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3.7 Varia-tion of 

^ o 

Shells viTith different geometric parameters are analysed 
by the pro-'Osed lower bound solutions* The values of are 

T> -p 

found for different ^ ratios with a fixed ratio epual to 

100 and for v.arious values of <f>^. The angle p which defines 

2 

the position of the neutral axis is taken as equal ^ 

The value of p ^ in the lower bound solution 3 is taken jS?om 

*R 

the Table 3.1 The variation of p^^ and v/ith ^ ratio for 
various values <1!^^ is plotted (Pig. 3.2), These graphs are 
useful to ostimate the ultimate load for any given sholl. 

3.8 Analysis of Shells 

T' 0 naells v/ith different geometric parameters are 

analysed by both elastic (26) and the proposed lov/er bound 

solutions, dlastic analysis is carried out for a wor]cing load 

of iJitencity 290 kg/m^. The ultimate load and the stress 

resultants ore calculated at collapse. In the limit analysis 

the amount ol transverse reinforcement is taken such that 

2 2 

11= 0.236. a end are equal to i50 kg/ cm md 2600 kg/cm 
c sy 

respectively. 

Per estimating the collapse load for a given thickness 

of shell, 'f v/ill be the total thickness of the shell while 

using K -yield criterion ond will be the effective thickness 
c 

while using cxiterion. 
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3,8.1 i-TuMerical Example i (Long shell) : 

Hie geometric parameters of the shell i are j 
P Tj 

Q = 30 , ^ ^ ^ 3 r333 , R =: 8 #0 HI aiid "b = 8 #0 CiU-S # 

Tliio snell fails due to I -yield criterion. Solu- 

c 

tion 2 gives tlie best lower bound and is equal to 640 kg/m^. 

Hence the loo-cl factor is equal to 2.2, The values of the 

stress resultants at critical sections at collapse are 

tabulated in Table 3.2, The stress resultants in reduced 

form at critica.1 sections are plotted both for elastic and 

limit analyscf, for a comparative study (Eig. 3,3). The 

IT H H M 

stress rcsulUnts m reduced form ere ^ and 

where == 0,11!^, The variation of the stress reriultsnts 

y. 

with arc [jlobted^ v/here y^ and li are the verfcj-cal distance 
measured Iron the springing level and depth of the shell 
respective I-' . 

3,8.2. Huiierical Example 2 (liedium shell) ; 

The geometric parameters of the shell 2 are j 

= 4 5°, p =1 <{>^ ^^ = 0.5818 , I = 2.5, R = 8.0 m end 
t = 8 ,0 cus . 

Tills shell fails due to yield criterion. The 

solution 2 gives the best lower bound vhiicli is equal to 
560 kg/is^ and hence the lx>ad factor is 1.9 5. The values of 




Table T ,2 

- Reduced 
sections 

stress resultants 
for shell 1 (Ion 

at critical 
,g Shell) 

<(> 

X 

% 

(x = L/2) 

4 * 

E 

0 

M 

0 

So 

(x = 0) 

0° 

-1.0000 

-0 ,.0697 

-0.7072 

0.0000 

5° 

-1.0000 

-0 .0670 

-0 .66 58 

1,0476 

o 

o 

-1.0000 

-0 .0 589 

-0.5i93 

2 .09 52 

15° 

-1.0000 

-0 .04: 56 

-0 .3808 

3 .1428 

o 

o 

0 .00 oc 

-0 .0271 

-0 .1985 

4 .1889 

2 5° 

1.057'-' 

-0 .0078 

-0 .0548 

3 .2082 

o 

o 

4.266,’ 

0 .0000 

0.0000 

0 .0000 



Table 1 .3 _ Redaced stress resultsnts at critical 
sections for sliell 2 (Hediuai Shell) 




<l> 

I'i 

(x =. 1/2) 

i'o 

3i 

J.U 

0 

K A 
x*!) 

S 

0 

(x = O) 

0° 

0.1127 

-0.0457 

-0 .9376 

0 .0000 

5° 

0 .1-63 

-0 .0469 

-0 .9340 

-0 .2936 

O 

O 

-0.1100 

-0.0496 

-0 .8900 

-0 .3329 

15° 

• 

c 

I 

-0 .0 5l4 

-0 .8094 

0.0477 

o 

o 

-0 .G(j89 

-0 .0492 

-0 .6849 

0 .8393 

25° 

-0 .Gove 

-0 .0408 

-0 . 5164 

1.7793 

30° 

0 .000 7 

-0 .02 72 

-0.3304 

2.2825 

3 5° 

0 .34 iG 

-0 .0131 

-0 .1596 

2.0487 

O 

o 

0 .728 i 

-0 .0028 

-0.04i6 

1.3064 

45° 

1.3000 

0 .0000 

0 .0000 

0.0000 


Tcblc 3,4 - lower bound solutions 


Shell 

Type 


’ L oT/er" iVund” 

3 oil t ion 
(p jn kg/m2) 

Lower bound 
solution 2 

(p in kg/m^) 

"l Ou'er* b ouiid 
solution 3 
(p in kg/m^) 

c 

Crite- 
r 1 on 

Crite- 

rion 

I 

c 

Crite- 

rion 

^ if) 

Crite- 

rion 

L 

c 

Crite- 

rion 

Crite- 

rion 

H ^ ^ 




L-**- • — -- 

* - 


lest lower 
b ound 

(p in kg/m^ 


Long 538,0 851 640 »0 056.0 540.0 960 .0 64o .0 

Shell 


Medium 
Shell 2090^0 


320,0 1382.0 


420.0 8l2,0 560.0 560.0 









FIG. 3. 4 DISTRIBUTION OF STRESS RESULTANTS-EXAMPLE 2 
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the stress resultants at the critical sections in reduced form 
are plotted for hotn elastic and limit analyses for purpose 
of comparison (fig. 3.4), The stress resultants are tabulated 
in Table 3 .3 . 

for a comparative study the values of the lo»/er bound 
computed from the three solutions, are presented in the 
Table 3.4, fox the two shells. In each solution the collapse 
load p is computed using both It and M -IT yield criteria 

C (J> <j) 

and the smaller of these tv/o values is taken as the loafer 
bound. In like manner the lower bounds are determuied by the 
three solutions and the best lower bound is indicated. 

3,9 Design of Shells 

In this section the design procedure for a simply 
supported DC cylindrical shell roof with free longitudinal 
edges IS described. Two shells are designed one for each 
mode of collapse. Elastic design is also carried out for 
w or k mg 1 o o-d s . 

Tho step by step procedure of (design is is follov'rs : 

(a) For a given vrorking load, the design ultimate load 
will bo estimated using suitable load factors. For 
simplicity of calculations, distribution of live load 
can bo taJren as gravity loading. 
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(Td) The value of p will he preassigned. value will be 
taken Cron Table 3.1, 

( c) The noiidinensional loo,d will be estima.ted froa the 
three lov'er bound solutions or from graphs ('n’li. 3.2). 
Out of chese three values the maximum value will be 
picl^ed up. The lower bound solution chat yieldn the 
maximum vclue of p^ will be noted. 

(d) Knowing p^, the ultimate load v^ill be calculated from 

the relation Pq = ^ where p is the ultimate load. 

c 

Depending on the relative values of the calculoted and 
design ultimate loads o and t v/ill be revised. 

(e) The stress resultants will be computed, using the 
corresponding lower bound solution, 

(f) The ti-onsverse steel \/ill be found from the Bq_, (2,29) 
so as to follow the force distribution at different 
points on the shell. 

(g) Longitudinal and diagonal reinforcement v/ill be computed 
from tne distributions of and respectively. 

In ulL these calculations, the stress in the steel 
will be tliG yield stress ^ 

To compare with the elastic design, iv;o shells with 

different geometric parameters are designed. Llastic design 

2 

IS carried out for a total \/orking load of 290 kg/m and 
limit design is done with an average load factor oC 2 ,o . 
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Shell 1 fails due to N -yield criterion and solution 2 gives 

w 

the "best lov/er hound. Shell 2, fails due to yield 

criterion a,iid solution 3 gives the best lov\/'er hound. The 
parameters of the two shells are ; 

Shell 1 (long Shell) : = 30°, § = | = 3.333 and 

E. =: 3^0 

Shell 2 (kodium Shell) i <()^ = 4 5°, ^ ^ p ^ = 0.38l8, 

~ = 2.5 and R = 8 .0 m . 

For ti'c tvJo shells, thiclcness of the shell, compressive 
strength of cencrete and the quantity of steel at critical 
sections a,ro presented m the Table 3.5. The qu 8 ntit 3 - of 
longitudinsil steel snovm indicates the steel required to 
resist the oo uOl longitudinal tension developed at the centre 
of the span. The transverse steel is estimated at crovm 
to resist xi^ end per unit length. Diagonal steel per 
unit length lo computed for resisting the diagonal tension 
equal to h,, x developed at support. 

T 

Prom the Table 3^5 it can be seen thot for shell l, 

v/hich fails due to N -yield criterion there is a considerable ' 

c 

economy due to the reduction in the thichness and strength 
of concrete, m the limit design as compared to the elastic 
design. There is a saving in steel of 15 per cent in the 
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Table 5 , 5 - Corfloarison betv/een elastic and limit designs 



Till cine SG 
t in fi'is 

in kg/ cm 

1 Transverse 
steel in 

2/ 

cm /m 

Longitu- 
dinal 
steel in 

2 ; 
cm 

Diag onal 
steel in 
2 

cm /m 

Shell 1 
(long Shell) 





a) Elastic 
de s ign 

10.0 

200.0 

8.8 

102.0 

11.5 

b) Limit 
design 

8.0 

150 .0 

6.0 

86 .0 

17.0 

Shell 2 
(lledium She’’ 1) 





a) Elastic 
de s ign 

10.0 

200 .0 

9.8 

- 51.0 

8.6 

b) Limit 
design 

9.0 

l50 *0 

11.0 

65.0 

11*5 
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longitudinal and 33 per cent in the transverse directions 
hut the increase in the diagonal reinforcement is 30 per cent. 

Tor the shell 2, 7hiich fails due to yield 

criterion, there is an economy in the concrete due to 
reduction m tne thndcness of shell and the strength of 
concrete over elcistic design. But in the limit design 
there is a conorderahle increase of the amount of steel. 

The increase is 27.5 per cent in the longitudinal direction, 
l2 per cent in the transverse direction and 3 5 per cent in the 
diagonal steel, Jlence the total economy depends on the 
relative costs of the materials. 

These figures gives only a comparative idea hut not 
absolute values as these are computed at the critica-1 sections 
only. 

Prom the foregoing analysis of shells it is concluded 

that long sholls fail due to the B -yield criterion and 

o 

short shells fail due to Jh-N, yield criterion. 

4 > 9 



CHAPTER 4 


LIMIT ARiLYSIS OP CILINDEICAL SHELL ROOPS 
\?ITH LORGITUDIIAL ELOE BEAMS 


4.1 G-eneral 

The load ca-rrying capacity of a shell can he increased 
either hy proviJing edge heams or hy increasing the thickness 
of the shell. In the case of short shells the latter proposi- 
tion is economical while the former proposition for medium 
and long shells. This is more so in the case of long shells. 
With this in vicv/ lower hound solutions are derived for 
simply supnorteci RC cylindrical shell roofs with longitudinal 
edge heams subjected to uniform gravity loading. 

Three lo'/er hound solutions are developed using the 
same procedure presented in the previous chapter. Tv/o shells 
one long and other medium, with different geometric parameters 
are considered. The elastic axial ysis is also carried out 
for w'orking loads. The distribution of stress resultants at 
critical sections are presented graphically for a comparative 
study. Two shells are also designed hy both elastic and 
proposed methods and relative economy is discussed. 

4.2 Lower Bound Solutions 

The lower hound solutions for the class of shells under 
consideration have to satisfy the eq_uilibrium equations, 
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Ijoundary conditions and the yield conditions. The equilibrium 
equations axid yiold conditions are presented in Chapter 2 in 
Sections 2,2 and 2,3.6 respe ctiyely . The stress resultants 
developed m the shell have to satisfy boundary conditions 
at the 3 ’anction of the shell v/ith the edge beams. In other 
v/ords the interc.cting forces developed at the interface of 
the shell and the edge beam should be compatible. 

Consider a RG cylindrical shell with edge beam. The 
coordinate system is shown in Pig, 4,ia. The geometry of the 
shell and loading are symmetric about the X axis i.e., at 
= 0 and the longitudinal edges are stiffened with edge 
beams at <? = Here is fcne angle measured from the 

cro'i/n to any ooint on the shell across the section and 
is the semi central angle, Tne cross section of the shell 
13 divided into t\io zones, viz. the too compression zone 
for 0 < < p and bottom tension zone for p < <P v/here p 

is an angle dcfimmig the location of the neutral axis. The 
stress resultaiitc have different expressioxis in xlie tv^o zones 
hut have to satisfy the continuity conditions at the region 
boundary of the t\ro zones. 

4.3 Equilibr luru of the Edge Beam 

The edge beam is in equilibrium under the action of 
its self weight, the longitudinal force de^veloped in it and 
the interactmig forces at the interface of the shell and the 
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FIG 4.1 (a) COORDINATE SYSTEM (WEDGE BEAM SHELL INTERACTING 
FORCES (c)EDGE BEAM FORCE DISTRIBUTION 
(d)EDGE BEAM LOADING 
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2H = depth Ox the edge heam, 

27/ = V7idth of the edge "beam, 

From Fi'_ , 4,ic it can he shown that 


a = 


^(2 ^n) 
5( 1 + n ) 


Talcing the vcnotion of F in the longitudinal direction as 


X 


parabolic. 


-S'. 


^ (x) (L _ x) 

.. L 


• • • (4 • l) 


At any distance :: from the origin, the equilibrium of the 
edge beam in the longitudmel direction gives that 


/ 


X 


IT , dx = 

ZZ9^ X 


... (4.2) 


Differentiating both sides with respect to x 


d F 




X 


dx ’ 

4 F^ 

_ —"FT” ( 1 — 2 x) . 


... (4.3) 


Let and be the resultant vertical and horizontal 
components of the interacting forces acting on the edge 
beam. They are given by 


.= 11 sin ^ + Q cos (? , 

V <^0 ° <{>0 ° 


V, =. IT. cos 0 - sin <P > 

h <l>o o 'J'o ° 


... (4.4) 


... (4.5) 
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Using the assLuaption (ii) Bq* (4*5) yields 

Q = U , cot 

Eqs. (4.4) and (4.6) give 

U = V sin <l> • 

0 

From. Eqs. (4.6) end (4.7) 

Q = cos . 

Now the edge beam is m equilibrium under the 
of a N and T (Eig. 4.id) . Taking moments about 
of the edge beam, external bending moment ll at x = ^ 

9 9 

q.1^ 

8 ~ 8 ’ 

and the internal resisting moment at x = ^ is 
= Eq a H . 

Equating Eqs. (4.9) and (4,lo) one gets 

8 E a H 

V = (q — -°' 2 ” ) • 

Substitution of Bq. (4.1l) in Lqs. (4.7) and 
give the forces IT^ and Q respectively. 


. (4.6) 

. (4.7) 

. (4.8) 

action 
the top 
is 

.. (4.9) 

.. (4.10) 

... (4.11) 
(4.8) 
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4*4 Symmetry, Boundary and Continuity Conditions 

I'lie shell IS symmetric m geometry and loading at out 
X-axis i.e.,at (*> = o» Hence only half part of the shell is 
considered using ciie symmetry conditions at <P = 0 • As the 
cross section of bhe shell is divided into the top compression 
zone and hottom tension zone, the stress resultants have 
different expressions in the two regions. But the stress 
resultants have to satisfy continuity conditions at bhe 
region boundary. Lastly, the stress resultants have to 
satisfy the boundary conditions at <!) = i.e., at the junction 

of the shell with the edge beam, fhey are summarised as 
follows : 


(i) Symmetry conditions : 


J'T 


x<l> “ 9 


9 9 


= 0 at <i> = 0 


» . • (4 « lE) 


(ii) Continuiuy conditions s 
I 


i'T, £j.icl ri, and the derivatives . 

X'l’’ 9 9 99 09 

continuous at 9 = S . 


M9 , 

and are 


... (4.13) 


(ill) Boundary conditions : 

The stress resultants and at 9 = 9q 

should be equal to H . , 11 , and zero respectively. 

^*^0 '^o '^o 

Therefore using Eqs. (4.3), (4.7) and (4*8) it can be 
shown tha.t. at 9 = 9 q 
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4 J 

= -^ (I - 2x) , 

**(, = t'o ’ (4.14) 

= \ OOS , 

where is giyen hy Eq. (4.ii). 

The shell is considered as simply supported at the 
ends. At the simply supported ends N = 0 at x = 0 and x = L. 

4,5 lower Bound Solution i 

In this solution the distribution for E is assumed 

as uniform across the cross section of the shell both in the 

tension and compressicon zones as shown in Eig. 4.2b. As 

E IS uniform across the section, it is independent of (ji • 

This kind of distribution of E can be expected for long shell 

where the behaviour is predominantly of the beam type. Yyhen 

the E -yield c3:itGrion governs the failure (collapse), then 
o 

the collapse mode of a lo/ig shell is similar to that of a beam. 

Sine© tne magnitude and longitudinal variation of E 
IS unkno^'/n 5 it is evaluatec' from the static equilibrium condi- 
tions of the shell, considered as a simply supported beam 
subjected to uniform gravity loading of intensity p. Then 

the other stress resultants N„ E , Q and M are obtained 

^9 9 9 9 

using the equilibrium equations (2.15) through (2.i8). 
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Tlie loii^ ibadixial force developed in tbe edge 163331 is 
taken es unifoi i across the depth. In other i-zords steel is 
distrihutel unifonily across the depth of the edge beam to 
produce uniforii force per unit depth. Hence q is taken as 
equal to i (fig. l.lc). 

4,5.1 Stress resultants 

(i) Evaluation of E j 

let 'Ejj^ and represent the magnitudes of m the 
compression and tension zones of the shell as shown in 
Eig. 4.2b. Inus 

for 0 < <J) < ^ , (4.15) 

I'T^ = for p < . ••• (4.l6) 


Since If. IS uniConn across the section, the resultant 

compressive anf tensile forces act through the centroids of 

the correspondia parts of the shell on either side of the 

neutral axis. I'St h. ard ho be the vertical distances of the 

1 ^ 

centroids of compressive and tensile regions from neutrei 
axis respectively. They are given by 

sin p 


h^ = R 


hg = R 


cos p 


cos p 


( sin - sin p ) 

IT 

^ ^0 



Addition of and lig 


gives 


-^2 = 

where 

C 6q sin p ~ ^ sin 

While consider mg the static equilibrium of the shell 
tensile and compressive forces are taken as positive and 
negative respectively. The static equilibrium condition 
(algebraic sum of the forces in the longitudina.1 direction 
IS equal to zero) results in 


2E 


XB 


li( - p ) + R p = 0 


( 4 . 17 ) 


where R, is the longitudinal force in the edge bean, 

Tne mouiexi-c: of resistance at any section x is obtained 
by taking monexits of the internal forces about the centroid 
of the compresGiOxi zone. Thus 


il 


r 







2f 

X 2 


R, 


where 


a 


2 - 



ii 


oos q + 5 


The extornol bending moment M at any section x subjected 


to uniforml^A distributed gravity loading of intensity p, is 
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The second tern on the right hand side is due to self weight 
of the edge heo is. 

Equating blie moment of resistance and the external 
bending marent II axid substituting E from Eq, (4.i), the 
following result is obtained 


N. 


XB = 2( ^ R *^0 i ” ®^o - x) . ... (4.i8) 


Substitution of Eqs , (4.l) and (4.18) m Eq. (4.1?) yield 


N, 


XA 


v-rhere 


2^a R 


(P 


0 ^ R 


L 


3^-0) 


... (A. 19) 


mt aawi’TJwrafw.wniiwiirii'inii iw 'iwni m 

3 ■ ( u - FI 


Jii 

;in _ sin p _ ( (j _ p)(cos cj) - ■^) 


0 R' 


Prom Substitubion of Eqs . (4.19) and (4 .i8) in Eqs . (4.i5) 
and (4.i6) respectively ib is obtained 




‘R ^o + 1 - ^^-3 

1 ij 

for 0 < <!> < P 


... (4.20) 


IT = 


X = 2TT“ - y)a R <^0 R 


8”^ cl 

^ (p S '* •^^“) (x) (1 - x) 

for ^ 4 '^’q • 


Evaluation of IT 




.. (4.21) 


The value of ^ is obtained as follows : 
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Eq. (2, 15) can "be written as 



R 


dR 

ix 


X 


, (4,22) 


( i) Eor 0 < ^ 


Substitution of Eq. (4.2o) in Eq. (4.22) yields 


dT^ = "^0 + R " (i - 2x) . (4.23) 

Integrating vvitli respect to one gets 

8cl !F 

“ S^aT" '^0 1 " (I - 2x) f D(x)j ... (4.24) 

^ '1 Ij 

v/here I)(x) is sn arbitrary function of x only. 

(ii) Eor p < <> < 

Substituting Eq. (4.2i) in Eq. (4.22) it is found 


= - Erfr-'^?')U t -p I - d- - zx), ... (4.25) 


which gives on integration ^-^ith respect to 


I 


X <!> “ 2^ $0 - pi £U, *^0 


(L _ 2x) + E(x), 


... (4.26) 


wnere B(x) is an arbitrary function of x only 
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Evaluation of li a ; 


In IS c'eteraiined in blae followin-g nanner 


]q. (2.i6) IS rearranged as 


2 


o 0 N 

v2 X 


IT. + ■ = E 2pE cos (ji 


d (J> 


dx 


. .. (4.27) 


(i) For 0 < «!) < 


Substituting for JM (Eq. (4.20)) 

Q 8a„ P 

IT^ + — --g™ = g'*’^ (p cJiq — 2pE cos <i> ♦ ... (4,28) 


d 


p a. 


'o ^ E 


The foregoing differential equation has the solution in the 
form 

p 8a„ P 

N = A sin i{> + B cos <;) - pE <> sin (J + '^’o R “ > 

V 1 -1- ° 1 B 

... (4 . 29 ) 


where A and are constants of integration. On differentia- 
1 1 

tion Eq. (4,29) gives 


I ) 

^ SB. 

= A, cos - B . am ^ - pR(sin + <> cos <>) . 
5 9 1 1 


... (4.30) 


( ii) Eor p 5. S. **’ o 

Eqs. (4,2i) and (4,2 7) give 

... (4.31) 
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The solution to the differential equation is 


E 

= Ag sin + Bg cos 4 > - pB 41 sin (ji - (P ***0 + 


4 -0 . 


... (4.3?') 


v/here Ag and Bg are constants of integration. Differentiating 
one gets 


■^ = Ao cos 41 - Bo sin 41 ~ pR (sin <> + 4 ) cos 4 )) . 


.. (4.33) 


Evaluation oC II||, : 


I'l is evaluated as follows : 
4> 

Eq. (?.1V) IS rewritten as 


E — pR COScj) 

d4>‘^ 


... (4.34) 


( 1 ) For 0 < ^ < p 


Suhstiaution of Eq. (4.29) in Eq. (4.34) yields 


5 ^,'i 2 

= -A R sin ~ B. R cos 4 . + pR ( 4 sin 4 * - cos 4> ) 

6 p 


Q 8 a^ F 

/ 4 o Os 

(p ^ - -^p—) . 


... (4.35) 


On integrating once it gives 
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6 M^ 

bT 


A. R cos <;> - R Sin <) _ pR^ (p cos P - r“'“(p + i 

1 1 ^ 8 it 

1 


8 a„ 

O 0 \ ^ 

— T-) + , 

Jj 


... ( 4 . 36 ) 


where is a constant of integration. Integrating once 
again 

2 

R sin 4i + R cos 6 - pR ( cos 6+4’ sin 6) 


^ 6 ^ . ^ ®^3 ^0 

2p'‘a~ 0 ■*' R ■* 


+ 0 6 + Cg , ... (4.37) 


where ^ constant of integration. 


( ii) For p < ‘^’ < 6^ 


From Eq.s . (4,32) and (4.34) it is found 


a*" 


2 

= -AgR3in6 -BgR cos 6 + pR (4 sin 4 - coo 4 ) 


. 4 lO) 


... (4,38) 


Integrating once 


56 


R^ 6 


= Ag R cos 4 - Bg R sm 4 - pH 4 cos 4 + a”"^P 6 


4 ®®-2 ^Ov . 

+ R - - T— “> S 


... (4.39) 


where is a constant of integration. Integrating once again 

O 
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= Ag R sin <) + Bg E cos (ji - pR ( cos 4> + <P sin ^ ) 

rjS ,2 Q Sap F 

+ (p (>^ + ^^ - + C3 <> + C^, . .. (4,4o) 

0 ” 1 Jj 


where is a conrtant of integration. 


Evaluation of J 


<!> * 


IS ohtained as explained below : 

<p 

Eq. (2.18) inay be written as 

Q. - 


I 

“ R 


• •• (4:«4l) 


Substitution of Eqs . (4,3-6) and (4.39) in Bq. (4.4i) yield 


R<j) Q 

= Ap COS <? - Sin 6 - pB C> cos $ - p~ (p f - 


3 0 




R 


-fov 


R<^ , ^ 

R 


Q<j = Ag cos <J> - Bg sin (> _ pR <> cos 6 + ‘^'o 

8ao E C„ 

- — ^ ’ -po/ (i £ t ^0 • • • 

L 

respe ctivelj'", 

fhe coxistants of integration are written in the 


modified form cs follows : 


Bo E 


A, = A^p pR+ A^g q + 8A^3 , 


‘12 


13 


B 


Bo R 


^ = B^^ pR + Bpg q 4- SB ^3 , 


Bq R 

Ag = Ag^ pE Agg q 4- SAgg 


5 
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E H 

Bp = Bp pR + Bgg q 8 3^3 ^ , 

^ ^ Jj p 

°1 = °11 * ®12 ® °13 ’ 

-o ^2 

? ^0 ^ 

^2 = ^21 *" ^22 ^ ^2S ’ 

^0 

C3 = C.^ pR^ + C32 qR + 8 C23 , 

2 ^0 

pR + C 42 qR + 8 0^3 . 

I'lie constonts of integration are determined from 
the conditions (4,i2) tlirouga (4.14). Thus 


I)(x) 

E(x) 


= 0 




iT 


(p ^0 ^ S 


8E Up 4B 

(L- 2 x) + (B 

1 B 


^11 "" "^12 " "‘1^ " ° 
^ s in B 
A = -° 

^21 a Q 


A, 

^22 = 4) 


21 


^23 " a 


0 

sin ^ 


( 4 ) cos <!>„ - sin 
^0 o 


qp ^ '*'0 

4 ) (p cos 4 q - <i>Q cos p ) 


■£'i- - 0 

11 '-^qf \ '<'0 


B-- = — -anTA'-'iT- 


B - Bj, 
®12 ' 


^of) 


2x) 
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®13 = - - <^0 ^0 ^*^0 - 


sin ( ^ ; 


B 


cos 

_ o o 

21 - S^<j~‘l"|7 


B 


22 - (j) 


21 


B 


23 


(% S cos <>^ + sin 6 COG <> + (i> sin 

ii ' 0 ’ O O 

P - p) 


sin 'J’q - P ) 


^11 = °12 = C^3 _ 0 


0, 


21 


«o (2 + P) 

cos 


'22 


(^P + 2) 

"#a7— ^1" ^ 


'23 - - p)p 


O ^ 

( sin p sin (}ic + 2 2^ 


H 

1 ^ "0 


+ ) cos <!)^ 


sin ❖ + sin 

0 


, ^2 >, P 

- <) ( 1 - 2a " 2 
1 


'31 = a:TTrrn 


'32 


“1'^ 0 

±0 

^ <J>o “ 


"-2 


'33 = ayT;i|i 


^ 1 
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°4i “ [ '^o =“ P “>0 - ^ ^ 

°42 = aS^pT^'^T L '^O ? - 8«(, ®“ P S3J1 (,^ - p(l + cos 2 «^) j 


°43 = S^f?;'^') [/“ 5 '>0 P 

,2 

P . n 9 a ^ 'H ^'2^0 

+ 2^1+ 20^))J - 

The stress resultants are summarised after substituting 
the constants of integration in the modified form as follows : 

( i) Por 0 <J> ^ 


'^x = 2pi^ <i> + I - "--2"-) (ri - =^) 


• « • *42) 


2pa^ 


a 

(p ^0 n ~ (L - 2x) 

1 


... (1.43) 


= -pE(« sill 0 -. p- -Bii * <l(Big =0= « + pi") 


8P Pu 


°2:'" ^^13 ^ 


... (4.44) 


p P 

- pR^ + cos <> + >{> sin ^ - G2-]_ ” ^11 


4® %-q - ®12 * - °2S^ 


P 

8R R^ <J>'' a 
0(6 


^13 ^23^ 


... (4.4 5) 
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= -pR cos (> 4 - sin <;>)-. q + B^p sin <?) 




^a. ^ ^12 


SB R a„^ 

(fa;' - ^13 


... (0-.46) 


(ii) For p 


R ®^o ^-2, 


(p « + 


I 


X 


2(0 -p) a R 


(x) (L - x) 


... (4.47) 


B 


x(J 


(P 0^ + 


q. ^2> 


o R 




2a p( 


(00 -«) + 


4B 


(n - 2x) 


. . .(4 .48) 


V 

= - pR ( 0 sin - ^21 ~ ^21 ^ 


+ q(J^22 * ®22 0* " ’fTTQ 


8Bo R 


a. 


•a;(T;;. pT 


+ Ao- 21x1 ^ + -Bg^ cos 0) ... (4,49) 


^2 , 

0 V. 


Ma = 


2 ‘ ‘ o 

-pR (cos « + <) sm <}> - 


- Ag^ sin ^ 


- Bg^ cos - C^^) + qR(A 20 sin ^ 

02 , 

+ Bgg cos 0 + + S 2 ^42^ 

8B R^ 

+ '“’§""(■^3 ^ ^ ®23 ^ ^ ^43 

L 


ag r 

^'^(1^4 y'> 

... (4. a) 
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Q* = - ♦ =°‘-' ♦ - ■ ■^21 ^ ®2l * " °3l’ 


+ dCAgg t-io « - Bgg sin ♦ + + O.g) 

1 0 ^ 

QjJl “PJ ^ 

+ (Ag^ cos <1. _ ^23 sin ^ + G33) 

... (4.5l) 

Prom tlie expressions for , II^ and Q^ it can be found that they 
are independent of x. 


4.5.2 Deternunrtion of Limit Load 


fhe limit load is evaluated by using the yield 
conditions given by Eqs . (2.26) and (2.50). Leu p^^ and p^g 
represent Lie o^'■o lower bounds obtained from the IT and 

w 

M,-iT yield criteria respectively, 'fhey are evaluated as 
9 4> 

follov'7s : 


( 1 ) Due to IT ^yiield criterion. The msximuia value of IT in 
compression 13 de ae mined frou Eq. (4.42) which occurs at 



IT 

X max 




8P^ a„ 
0 3 





8 p a^ R 


... (4.52) 


rrom the yield condition ( 2 . 26 ), it is found 


P 


11 ' 


8a p RIT^ 


^ L 
o 


1 


o 3 

<f>0 L 


... (4.53) 
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(ii) Due to Djj) - yield criterion : llie failure of a 
cylindrical shell occurs in tie nrensverse direction, if the 
yield condition. ^ lyen hy Eq. (2-30) is satisfied. In a simply 
supported cyliiuhrical shell ’ a.th edge ‘beau.s, li^j, is alv/ays 
compressive our may change in sign vritli in the range from 
0 xo There/ ore the yield condition is xo he cneohed v;ith. 

in the region 0 co Since and are independent of x, 

it can he understood that if once yielding is initiated at a 
point, failure occurs for the entire length of the shell. In 
this limit analysis of shells it is observed tnat when the 
shell fails due to the yield criterion, the failure 

of the shell occurs at the crown. The absolute maximum value 
of li(j, also occurs at cro^m. Consequently the limit load due 
to this criterion can be found using the values of and 
at crovm. let and represent the values of and 
at crO'wn, which are obtained from the Eqs. (4,44) and (4.45). 
Thus 


F 


<{>c 


= pF(; 


+ B 


11 




8E R 

+ 2— • (B 

1^ 


l3 ^ a. 




(4.54) 


M 


<J> c 


= pR + *^21 ■" i) + Gp2) 


11 


12 


'22' 


8F^R^ 

+ + Cgg) 

jj 


Substituxion of F . and M m 

<po $ c 


... (4.55) 
Eq • (2 .30 ) 


yields the collapse load 
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"12 


i- 

^1 


, '1®- ^ N cq 

b ~ r.— (B ‘ ^ ^ 


12 ^ ^22) fr + \2> 


LI^L' 


2 




3c P^ R 


V ^®i3 - ea- 
n_L r^i. 


^ry »««t 
® •) 


O 


... (4.56) 


wn re 


D. 


S: (P + c 

IIq '^^ll ^ ^21 


s cH ^ 
1) - fT" (S 
o 


+ ZO_S 

Ij - "“ll 


The collapse load p xs the smaller of the t’vo yalues 
p^l and pj^g. 

4.6 Lov/er Bound Solution 2 


in c.:o lower hound solution 2, the distribution for 
I IS as siiov/n in Pig, 4.2c, In this, is taken as uniform 
across the cross section in tiie compression zone and it varies 
linearly with oepth a' from the neutral axis in the tension 
zone. The reasuxis for selecting this distribution are same 
as that lor choHo v/ith free lont', itudinal edoes (Section 3.5). 

Tne venation of force developed in the edge beam 
IS trapezoidal as indicated in the Pig. 4.ic» The magnitude 
of the force per unit deptn at top of the edge beam is n times 
that of the one at the bottom. 

4.6.1 Stress resultants : 

Tne stress resultants and are 

determined as follows : 
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Evaluation of 11 ; 

X 

file iiia^mtude anci longitudinal variation of is 

Ia. 

found from, the occtic equilibrium conditions of bhe shell. 

let denote tno vertical distance of the centre of 
total coiipressio.i from the neutral axis and it coincides with 
the centroid of the circular arc of the shell above neutral 
axis. Thus 



The dep fells of the compression and tension zones are 
represented by and hg respectively. They are given by 

h ^ = 1( 1 - cos p ) 

hg = E( cos p - cos . 

If a* ic the -Terticol distance from the neutral axis 
to any point on fclie cross section of tne shell, then z’ is 
given (Fig. 4, 2a) as 

z’ = Pl( cos - cos p) for o <. * £ 

z' = R( cos p - cos Ip) for ^ ^ 

Tne I'lO-mmuiii value of 11^ in the tension zone and the 
uniform value in the compression zone are represenied by 
and respectively (Fig. 4, 2c). Then the distribution of 
in the two zones is expressed as 
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z' IT^(cos^ - cos <j)) 
2 ( cos p - cos <{>^) 


for 0 < 4' < p , 

... (4.57) 

for ^ < <^ < 'l>Q 


The force in the edge 56321 P at any section x is given 

X. 

5y Eq. (4.l) . 


Talcing tension and compression as positive and negative 
respectively, the siatic equilibrium condition of the shell, 
i.e., L E ds = 0 gives 


2 f E, 


X 


d <?) + 


I 


R d<J) +• 2P = 0 * 


Substitution Oj 


II 


XA 


R 


Sqs. (4.57) and integration yields 

(bsP> -1 

(4>o - <^0 + P J 


( cos p - cos ^Iq) 


+ E.. = 0 


from which the lollowing is obtained 


II 

P 


'^Sla f 

— ' ’ ( cos B 

2 


cos 


t) 


T 

_x^ 

Rao 


( COs p 


cos 4 q) • 

... (4.58) 


Y;here 

ap = (4^ ~ p) cos p - sin + sm p . 

Talcing 11 orients about tiie neutral axis, the moment of 
resistance is obtained* Tuus 
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^ o 

LI = 2 j IT 11 z’ + 8 / TNT R d4> z' + 2P^ (h„ + sB) . 
o P X 

Substidutiori. of from Eqs» (4»57) in the foregoing equation, 
gives 


IT 


r'- 


iA ® 


+ 2K 


X 


where 

= (sin ^ - p 
+ ^ (sin 2 


cos p) + 




~ sin 2p) ^ 


ag H 
(♦o-f ) 

^ — (2+ COS 2p) 

2 cos p(sin - sin p) 


(4.59) 




= ( cos p 


cos 


<^0 + 





(2 + cos 2p) 


+ ^ (sin 2 - sin 2p) - 2 cos p(sin (j)^ - sin p) 

The external Toending moment IT at any section x is 

IT = (pR 4 )^ + q) (x) (L - x) . ... (4.6o) 

Equating Eqs . (a, 59) and (4.6o) and substitution of P.^ from 
Eq. (4.1) ono gets 


X4 


^X. 
2a X 

1 


(P 


1 

+ 1- 


0 0 


) (x) (L - x) . ... (4.61) 


The first pc^rt of Eq. (-1.57) and Jiiq. (4. 61) give 

Q Q’F B^rr 

^'^x = - Ta"^ (P ‘f'o R “ (x)(l - x) for 0 < ^ < p . 

... (4.62) 



Il6 


Similarly usin^ Eqs . (i-.l), the second pert of Eq. (4.57), 
(4.58) and (4.6i) it is found rhat 

0 q 81 

R - P 

1 ^ L 

for p < ‘i’ < <!>q ... (4 .63) 

\;here 

H - f~ + ^-s- 

Evaluation of : 

IS determined as fo31ov/s : 

Eq. (2.i5) yields 



6 <() 


R 


dx * 


. .. (4.64) 


(i) Eor 0 < '(’ 0 

Substitution of Eq. (4,62) in Eq . (4,64) gives 




(p *0 + 


q. 

R 


810^3 

”Y") 
1 - 


(1 - 2x). 


On integration >rifch respect lo (j) 

n SR^a^ 

^‘x^ “ 2^ '^’o S *** 


where E(x) lo a function of x only 



llV 


( ii) Jor g < <1 < 4i 


?ron il'e Eq&. (4.63) and (4,64) it xs found that 


6H 


Xi(> 


5 <!> 


<^P *0 * E - 


- cos ) (Ii - 2x) 


On integration 


x<t> 


^ ^ A 

2a,a7 ^0 ■*■ R 

1 2 


8J'o^4 


^'^)( ^ COS p - sin<J») (L - 25 c)+E(x) 


• • . (4 • 66 ) 


vi/here S(x) is a function of x only. 


Evaluation of : 


E4. (2.16) iia^/ be rav'ritten as 

? 2 
5 IT. 2 6 N 

TJ, 4. — E « • ” '2*" — 2pR cos ^ • 

5 (j) 


... (4.67) 


( 1 ) For 0 5. ‘J’ 1 P 

SubstiGUtion of Bq, (4.62) in the foregoing eq^uation 

yields 


6\ 


U . “D 

= t: Cp '^0 1 


q ^’z 

• •^'— ) _ 2pR cos <}> * ... (4.68) 


+ ^ “ 


5(j)' ""I “ 1 

This differential equation has a solution ixi the follov/mg 
form 


sm^ + cos + I'’ (p '*’0 + J ~ — ^g-^)-pR <t> sin , 


4 ®^’o ^3^ 


... (4,69) 
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where and are constants of integration. On different iaticai 
one gets 


61^ 

sin ^ sin i)> ~ pR ( 4) cos <}> + sin <{i ) . ... (4.70) 

(ii) Bor ^ 


Substituting Eq. (4.63) in Eq. (4.67) 

0 31 a SB a. 

^4> ■*' "~2 = ~ T’aT" ip ’^o ■*' R ---^)(cos ^ - cos <J>)-.2pR cos 

6<{’ 1 2 “■ L 

• (4.7l) 


The solution oX the differential equation is 


= Ag sui ♦ + cos _ J™(p ♦„ + g - -^-)(oOS p 

1 6 ll 


4 sin 4 


) - pR 4) sin 4 > , 


... (4.72) 


where Ag and Jg ore coiistonts of integration. Different ia,t mg 
with respect to 4 > tne following is obtained i 


= -4 ♦ - Bg sin + g~a-(p ^0 + E - "75"’ 


p R 

T ^2 


8B a. 


( <J> cos ^ + sin O “ P^ ( ^ c; 0 S 4 + sin 4) ) . ... (4 ,73) 


Evaluation of s 

I'R IS evaluated as in tue following manner ; 
9 

Eq. (2.l7) gives 
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d 




o4>‘ 


~ R - pR cos <i> . 


(4.74) 


(i) Ror o < <J> < 


Substixutxon of Eq. ( 4 . 69 ) an Eq- ( 4 * 74 ) yields 




6<i> 


^ R 


COS Vq 


(P <>0 B 


Q. 0 3, 


+ pR <> sin cos $ , 


(4.75) 


Integrating t\/ice 




2 C ^Bq 3-3 

64) - -^1 B sin ~ ™ (p 4>o + R “ 

2 A 

- pR COS 4) + ^ • • • 


(4.76) 


and 


1 


2 2 

(, = H sui ^ + B E 00s « _ (p «„ + I -- qr ^ 


q BB ^ 


2a. 


- pR (cos 4+4) sin 4 ) + 0^4 + Gg, 

nrier® G^ ^nd Cg are constants of integration • 

(li) Tor ^ < <!’ < <!)g 

Prom Eqs. (4.72) and (4.74) it follows iBat 


(4.T?) 


^ 3 RT? a 

2 “ = —•^2 B sin 4 - Bg R cos (j) + ^"a" "' a' " o4^ 

(2 cOs ^ .. 4 31114 ) + pR^ 4 sin 4 - pB^ * 
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On integ rati on- 


ail 

^ 

5<!> 




q 


Ag R cos - Bg R sin (J + I " “TT^^ 

1 ^ 1 


(2(i)cos p - sin4> Hh 0 cos <> ) -- pli <> cos + 0^ 


... (4.79) 


where is o constoiit of integration. Integrating once again 
it is found tlio-t 

Ci Si’n^'4 

li^ = Ag R san(^ + Bg R cos « + (p *>0 R - ''7T~‘'^ 

1 JJ 

2 2 
( 0 cos §. + 2 cos -i- <) sin <>) - pR (cos <> + (Ji sm O 


+ 0 


3 (J 


4 - C 


. . . (4.80) 


where 0^ is a constant of integration. 


Svaluation of 'iij) : 


IS determined as follows : 
Eq. (2.l8) may he written as 



. .. (4.81) 


Substitution of Bqs. (4.76) a.no (4,79) respectively'' m 
Eq, (4.81) gives 


= A^ cos <) - B^ sin (t) 



61? Q 

^ (p g*") — pR 4' COS 6 

for 0 < <> < ^ , ... (4.82) 
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fi ?. , 

= Ag cos ^ - Eg sin ^ (p '}>o + R " ”^2 -* 


(2<^^cos p - sin ^ + <}> cos (j) ) - pR <»> cos ^ 


for p < <1’ < 


... (4.83) 


The const 


cuts of integration in the modified form s-re 


ore sente d as 


A^ = A^^ pR + + 8A^3 


Bi=BiiPR+Bi2 n_+ SB ^3^ 


\ " ^21 ^^22 ®^23 


Bg -- Bg^ pR + B 22 B + 8^23 p” 


Cl = C^^ pR^ t C^o oB + 80^3 


Og = Gg^ pR^'^ + Cgg pR + SCgg 


C3 = pR^ + C32 CB + 8G33 -p 

ToR^ 

O4 = G^^ pR^ + C42 + ^^43 • 

The consto-nts of integration are obtained irom the 

. 9 ^ +nrnneh f4.i4). 'Bisy given as : 
conditi-^ns (-±»l2) uui-ougn 
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D(x) 

B(x) 


= 0 


4? 


(P <^0 + R 


cos ^ 

Jj 


sin 


4- (L 2x) 


'^ll " " ^l3 " 


f„ + -- sin fi 

■^21 '^®-i^2 ^ ^1- 


^ In 
22 <^0 


-a, 


a_ 


A _ -* sin ^ 

■^23 " 4a^ag 3 a. 


3 




(i- 


li " 1 

B 


oos 2^ . 5} - ^ fg * *0 *< 


B 


11 


12 “ 4-, 


B 


^4^ 


13 


4a ^a^ *••"1 


(x - cos 2^ - 3j + ■^ -^2 


_ 3 ) + ^ ^2 ” ^ I 


1 


P' 


B 


21 " 4a ap ^l ■" a 


tan 4^ sin p + 4-^ tan <l>^ 


B 


B 


'12 


21 


22 " 4, 


o 


'/I P ^"3 

if -4- s in p 


Bgs = - 4a. So 1 a 


tan 't’Q - a ^ tan 


1 2 


Gii = ^12 = ^13 " ° 





Og, = =0® *0 + 2a 


^22 - - - - o 

^4 ^ j? -f -u a ^ &j-a <!> 

023 = - 4 + 2r ^9 * " E ^ 

1 

***0 ^ ^0 
^31 ^ ~ ^1 

. hi 


"32 " 

^0 



3*4 P 


^33 = 

” ^-1^2 

"^-l 


cos 4 'q 

Jsl 

11 

tM 

o 

+ 2a^a2 

n -* 

lol 

<H 

CO. i 

- sin 4> 

^42 " 

2 a ^ag 



a. ^ 


- -r- ^n ^ ^ ^ ^0 


<^43 - 2a^a2 ^6 ■" a^ ^o 
where f^, 3 - = 1,9 ere feiven hy 
f = tan <!>, C"- C. 0 S ^ cosec - sm 2^ 


_ 2 (4-, - ^)1 


f_ = cos ^ + sna ^ tan 

fCf 


sin 2^-2^ 

= sin ^ - P' eos p 

9 n ?* sin 4> ~ 2 cos <^ - 2 cos p 

cos p *" 2 <t>Q er Vq o 

+ 2 cos <i>Q 


^3 " 


fp; = 


= cos p (♦ o 


2 _ 2 * ^ + 2) + 2*0 
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00= g + (♦„ - ^ ^ *0 * 

_2p (cos + ♦„ "'o') - Z s“ "o + ® ^ 

f,= 2 (=osg - coo oosp -2(*„-^) 

o 

fg = - 2<^0 ^ ^ * 

The stress resultants are summarised arter suljstituting 
the constants or integration in the modified form as follows : 


( i) I'd 0 < 'S’ < ^ 


■fir Ti 

i‘x “ -a^R ^o R 


a (x) (1 - x) 


... (4.84) 


, 1 0 5 \ /y 2x^ 

_ _±- (n 4) + w - —“*^2 — ' ~ ^ 

^''x* " 2a^ ^0 R X. 


... (4.85) 


n. = pK Cr "r p°® ♦ - ♦ hs 


a. 


2 

9 . *^0 ^. . osr,R A _ 6 sin 4) + ^o'l) 

II = pR^ (B^^ cos - cos <? 21 


... (4.86) 


.2 „ ^2 


, .40,3) 

+ qR (B.e °°° ^ + ^22 ' 2a ^ ^ j^2 2a^ 10 

... (4.87) 


4)^4> 

-q/'d sin . 4 ■■■_'• 4 4* co*^ 

^ 

SR^R a 4) 

0 . - B „ sin 4>) 

5^ 1 fl. _ 10 


0) - sin ♦ + -;) 


1 


... (4.88) 


1 1 



125 


(li) For ^ ^ 


4 


®x = EiTi^E ‘^0 * 1 - -T~'^ ^ *) 

1 ^ 


• » * ^4 *89 ) 




a 8-^^ 3 -a 

H 0 4\ f“ 


= zi^ ^p '»o + R - L i.%-*)00S ^ - sin 


L 
4F 

+ sin<{>3 2x) + (1 - 2x) 

1 


... (4.90) 


^ Q P' ^ s in 

11^ = pE- 11^21 '*’ ®2i <>-'<*’ sxn i) - — — ' (cos — g~ — — )3 

1 2 


4’ sin 4' 


+ aC V SJE * + Bgg cos « - ^ (cos p - g )n 


8^0^ r- ^4 ^ / 4 S'3J1 4 \ -1 

+ “ir L-^2G + ®23 ^ ^ - ~S ) J 

L 1 


» • • (4 .9 1) 


M. 


sin 4 + coo O + <> :- - cos 4 - 4 sui 4 


21 


31 


^41 


^ B.(^) 

“Ta'a ^■^22 ^22 '*' ^32 ^ ^ 

1 2 


42 


8j? H 

+ ’-^1^) + (Agg Sin 4 + Bgg cos 4 + 4 + G^^ 

1 2 L 




2a^ag 


,«» (41: *92) 
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0 ^ 

% = ♦ - ^21 ♦ + Si - ♦ COS’ ♦ + srir^-’ 


? FgCO 

+ q.(Agg CDS ♦ - Bgj S3J1 * + 032 + 


83? E f, a ^ P2(+) 

+ -TT- (A03 cos ♦ - Bgg sm ♦ + O33 g rX— ^ 

J- 1 iC 

... (4.93) 


where 

2 

= <fi cos p + 2 cos <j> + <i) sm 4' 
P 2 (i{i) = 2 <|) cos ^ - sm 4>+ <> cos (j» 


?rom tno expressions for IT. , T.I. and Q . I't can be found 

9 4> <P 

that they are independent of x. 

4.6.2 Deternmation of limit load 


Tiie limit load p is the smaller of the two values 

and Pj ^2 e ''’‘o-loo ^e^l from the t >>'0 mdependent pield conditions, 

VIZ., IT and !I - yiold criteria re spe ctively. These are 
c 9 <p 

determined as follovie : 


( 1 ) Due to IT —yield criterion : T'le mci'-imum 

comp re sc ion v^rhich occurs at centre of bhe scan 

toll . It is deterr.ined from So. (4,84). Piu 
c 

1 ! ^ -.2 

U = — (p ij) + ft— “j 1] . 

X iiax 8a^E •^o E. -^2 
The yield condition (2.26) and Eq. (4.94) give 


voilue of IT,, in 

mi. 

IS limited 


. . • ( 4 . 94 ) 
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Pll 


# L‘ 
o 


-iL_ 

$0 B 


3?^ 

0 o 

. 

9 L 
o 


( 4 . 95 ) 


(xi) Di' 9 xo -q; - ''<u yi'-'ic’’ cxxteii .n ; As e:cpl:ined eorlxer in 
the Section 1.5.' i-ic fsilurc cf 5ae s^iell occurs in the 
trans-verse dirccxior ab tne croi/ii. Representing tne veiues of 
11^ and BIjj, at the crovn hy R and I'i (they are evaluated 
using Bqs. (4.86) and (4.87) respectively) one gets 




c 


M 


<(> c 


(j) 8R R a„ 

pR(^ + 3 ) + q(^ + Bq2) + “TT" *•' C4.96) 

1 - 1 1 1 

gy 

pR^(fll + - 1) + qR (B^g + Cgg) + (3^3 + Cg 3 ) . 

... ( 4 . 97 ) 


Suhstitution of Lqs. (4.96) and (4.97) in the yield 
condition giveji hy Bq. (2.30) give 


Pl2 = Up + c (-12 - lI7 (®12 + ^22^ 


0 


0 


8c R 


o 

+ “ "2“ 
H L 


^''l3 ” ^23^ 3 » *** 

o~ ^ ' o^ 


where 


R^ \ /T, On 

R = ^ (J + Cg - l) ~ (B + ^ ) 

0 o 1 


I’lie lov^er hound loa,d p is the spialler of xhe t'ffo values 
Pri ^32* 
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4,7 L ov/e r 8 ound S olut i on 3 

In tnc lo'/er Id ound solutions i and 2, the distribution 

of in comoression zone is token as nniform across the cross 

section. In hio lo^er hound solution 3, the distribution of 

N in compression is taleen as parch olic vdiich vanes with the 

vertical distciico z’ from the neutral axis. The maximum value 

of N in compression in the transverse direction occurs ot 

an angle measured from the corwn. The va^lue of p ^ is kept 

as a variable for a shell of given geometry and its value is 

evaluated such that it gives an optimal value for the ultimate 

load, Tne variation of H in the tension zone is kept as 

linear, same as 3-ii the lower bound solution 2, \;hich varies 

with the vertical distance z' beio\< the neutral axis. Tins 

€ 

IS shown in fit . 4 .2d . 

The distribution of longitudinal tension m the edge 
beam is taken bs trapezoidal as indicated m fig. 4,1c. 

The magnitude o* the force per unic depth at top of the edge 
beam is n times ti at of The one ax the bottom. 

The reasons for choosing ins kind of distribution for 
h are explamie. earlier in Section 3.6, wnile developing 
lower bound solutions for shells uith free longitudinal edges. 

4.7,1. Stress resultants 
Bvar uation of H., ; 

I'J is evaluated as follows : 

X 



1^9 


1 and II reD-’C'sent t'ne raaximum -values of N in the 

q p “ X 

compression an', tension zon's respectively as shov/n in 

Fig. 4. 2d. T'.e poi 'holic venation in the compression zone 

and the linesr variation in the tension zone are expressed as 


I'l 2’(2z' - z’) for 0 < <J> < S , ... (4.99) 

r7 t ^ 1 

and 

H z’ 

...(4.100) 

2 


where is the vertical distance from the neutral axis to 
the point of maicinium compression and hg is the depth of the 
compression zoxie. I'nus (Fig. 4.2a) 


z* = Ii( cos 'P - cos ^ ) for 0 < <J> < ^ ? 

z‘ = R( cos - cos *!* ) for p <. ‘J’ <, 

= E( cos ^ ^ - cos ^) , 
hg = R( cos p - cos <{>q) , 

Substituting these values in Bps. (4,99) and (4.100), ip is 
expressed as 

IT 

Pf _ _ »_iL . . cos P ^ - cos ^ - cos P ) ( cos P -cos p ) 

^ ( cos - cos p ) 

for 0 < <{' < p . • • (4.101) 

and 

K ( cos B - COG P ) 

for p £ 'i' £ 'J’o • 

^ (cos ^ - cos P ) 




(4.102) 
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virhere 


a. = 


C (<!>o ~ cos ^ - sui <) + sin p ;] 


1 r" 


[[2 cos p (sin p - p cos p) + 2 ^ U 


The moraent ol resistance at any section is found 
hy taking mocicnts about the neutrol aucis. Thus 


1=2/ NR a <!> 2 ’ + 2 / NR d4> 2 ’ + 2 R [/ R( cos p 


- cos + a H"1 , 
o 


... (4.108) 


After substituting Eqs. (4,10l), (4.102) and (4.10?) 


2 N R 


M = + 2^3 R , 

( COS p -cos p ) 


• (4.109) 


where 


rz Y 

ag = cos ^ . Up (2 + cos 2^) - ^ sxn 2p [] + J 2 ^ 


sxn ^ ^ ^ cos ^ cos 2^ + 


-f sill 


2^-3 sxn ❖ cos p + 3 sxn 2^]] , 


- cos p cos ^0 ^ R 


a _ Is. f 

R zi. L 


- S Sinci^o +“5 StriX^] 

The external bending moment U a,t any section x is 


r = (pr + q) (x) (L - x) . 


. . . (4.110) 
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Eqaatmg Eos. (4.109) axs-d (4, no) and substituting in 
Eq[. (4.i) for P„ one gets 


= gfcs Cp *0 + I - (cos - cos ^)^x)(L - x) 


... ( 4 • 111 ) 


Using liqs. (4.101), (4.102), (4.107) and (4,iii) 11, 


X 


IS expressed as 


®x = - (P « 


4 SS'oO'S, 


^ ^ ^ — )(2 cos - cos ^ ^ cos <!>) 

L 


(cos (p - cOs (x) (L - x) for 0 < ^ (4.ll2) 


H 


1 

1 


(P '^o + s 


^"0 

0 0\ 

~T““) 


2 a 


~c'r~ 


( cos s - cos 0 )(x)(Ij-x) 


for < <p <<i) , 


. . (4.113) 


where 


a. 


1 


n (<> - 55) cos p - sin e> + sin ^ ” 


Evaluation of 


:<P * 


U , IS determined as follows 

X(J) 

Eq. (2.15) gives 


x< 


6 P 


- R 


6U 

dx 


X 


.. (4.114) 
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(i) Por 0 < < {5 


Su'bstit'j.biOii of £q. (4.112) m Eq, ( 4 ,qi 4 ) yields 




d 2a^ 


q Sa^r 

'^0 R (2 cos _ cos ^ _ cos <;>) 


( cOs 0 - cos (ii - 2x) , 


... (4.115) 


winch on integration with respect to <{> results in 
1 <3 ®^3^o 

“xi. = Hr 4 + ii - " ♦ eos p) 


A Sin 2^ , 

I cos 2^ - - — ^ — ) (L - 2x) + D(x) 


... (4 • ll6) 


where D(x) is an arbitrary function of x only. 

(ii) Por ^ < 'J’q J 

rromfqs. (4.ii3) and (4.ii4) it is obtained 


JA - 


(p 

’®'1®'2 


4:3, S' 

(cOs p - cos <J>)(L-2x). 
L 


... (4.117) 


On integration \7ith respect to it yields 


^‘o, 


(p ^0 S “* "7^—) 

Jj 






( <P cos p-sin <> )(L - 2x) 


+ E(x) 


... (4.118) 
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v\'here E(x) xo c<si cr'bxtrsry function of x only, 


EvalUccXxo-i oi 




IS dete_ uned as described "below s 
Eo . (2.16) u re,7ritten as 


2, 

5 IE 


8 F 


2 '' r 

Na + — = Ji — EpR cos • 

d b 6x 


... (4,119) 


(1) For 0 < '^ < P 

Using Eq. (4.112) for 11^, Eq. (4.119) results in 


d^lT 

<>0 i - cos - cos p - cos b) 

0 c|) 2 L 


1 


(cos b - cos p) - 2pE cos b 


... (4,l2o) 


The differentxal eq_ustion has the solucion of the form 


sin b + E, cos b + *“> (p b^ + 
- 2 . os p 


^ s^AA^ Y 'r " a<^ 0 R 


1 


^”2) (b sin b cos p ^ 

cos 2s cos 2b, ^ , 

+ + — ^^) _ pR b sin b , 

... (4.l2l) 


ooc + -ng- 


where ji.^ and cro constonts of integration. On differentia- 
tion the above eqoo/Gion gives 


d b 


p 1 

^ A, cos b - E^ sin b + “ (p bq ^ 1 


Sa^P 

0 On 


1 ■ ^ ^-2 - p ' 

s 1 2 { 

(cos p^(sm b + b cos b) - ; - pR(sin b + b cos b)» 

... (4.122) 
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(xi) I'or p < 5 4 * 


Sub st ITU Lx on of Sq,. (n-.li 3 ) in Eq, ( 4 . 119 ) yields 


+ 


6 II, 


6 <}> 


= - R 


4 


(p 'J'o R ~ 4 a, R 


4^0 


®- 1^2 


( cos p cOs 4> ) 


2 pR cos (j> 


... ( 4 . 123 ) 


Tne solution of tne differential equation is 


IR 


Ag sm Ip + I >2 cos 


R 


Cp ♦„ + I- 


8 a 


3"o^ 


4 a, r 
4 0 


4> cm <f> 

( cos ^ - pR sin ^ , 


... ( 4 . 124 ) 


where Ag and Eg ore constants of integration. On differentiation 


-- — = A,, cos 4 > 
o ?|> ^ 


sin <^ + 


R 


q 8a,vR 

(P *0 i - --T-) 




( <J> cos +• 3111 <f>) - pR ( <!> cos <})4-Sin<|)). ... ( 4 ,i 2 5 ) 


Ewa-luation of li^ 5 


M. IS evalucted as follows s 

<p 


rom Eq. ( 2 , 1 ?) 


^ 2 

= ... bl R - piR COS 4> . 


or 


... ( 4 . 126 ) 
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(i) Tor 0 < ^ < p 


5qs. (4.1^1) dnd (4.i26) yield 


2 

piSr 

<P 

= 

5 


-A^ R sir. <;) — 3.R cos ^ 
1 ^ 1 


a ®Vo, 


(p ^ ^ 


( 4) sm (J) COG ^ ^ - 2 cos 


cos 2a cos 2<|) 

cos ^ ^ + — .g— + ) 


+ pR ( <!) sin <1) cos 4) ) 


(4.127) 


Integrating once 


t .2 q Sa^R 

A^R cos - B^R sin — (p <^o + ^ - 

o JJ 

(cos ^^(&i.i 4i - cos $) - 2t> cos cos ^ f 


sin 2<;) 


) -. pR <!> cos <() + C. 


... (4,128) 


ffhere 0^ is a constant of integration. Integrating once again 


M = A., R sin 4 - B R cos $ + (p + =d 
<{> 1 1 *^2 ^ 


Q 83 ,r 7 !P 

i 3 0^ 


(cos cos <> + ‘t* sm <!>)+<!> cos cos ^ 

cos 2^^ 2 

- 1“ cos 2^ + -■ ir“) - pR (cos (> + <!) sin <> ) + G^<J> 


+ Cq 3 


... (4.129) 


wiiere Cg is a constant of integration. 
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(ii) Eor p < <{) < (»), 


6% 

= 

6<J) 


Sabstifcutxng Eq_. (4 *124) xn. Eq[, (4.i26) 

(p + I - 


- AgR sin ~ Bg R cos <> + R 


<1 


^■1^2 


4ay,R 


p 

( cos e - ^ sin <^ ) + pR ( <!> sin tj) - cos (J ) , 

... (4.130) 


Integrating t^/ice 


dii 


$ _ X t.2 

: AgR cos !i> ^ BgR sin ^ + R 


Q 8a„F 

/ JL X ^ O 0 \ 

(p «’n + I - ~T-> 


^l®-2 


IR" 


1 ^ \ 
( ^) cos ^ ^ sm <) + “ cos « ) 


pR^<> 


cos <) + Gg , 


... (4.131) 


Ra = AoR oin <> t- BoR cos c> + R 

9 c 


, I 3 0\ 

<^0 "■ R ’• 7s -> 4a^R^ 


^1^-2 




( QOS ^ + cos *5* + 2 sin <> ) - pR^ ( cos <? + 3 sin <>) 


+ Cj (!> + f 


... (4.132) 


where C3 and are constants of integration< 
Evaluation of Q s 


IS determined as described s 
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Eq_. ( 2 . 10 ) oxy "be v^ritten as 


1 


- R sT * 


... ( 4 . 133 ) 


Substitution o:; 


-.IS. (4.128) ai>l (4.l3l) in Eq. (4.155) giTes 




00. « - 3, sxn , - (. «o . i - pqso. ♦ 


Sin 


2<>. 


- <;> COB O - 20 cos cos ^ + I cos 4 + -^2 ^ 


- pR 0 cos 0 + “ > 


... (4 *l 34 ) 


and 


Q = A„ cos 0 -Sg sin 0 + B. 
''A ^ 




(0 cos ^ ^ sin 0 + I cos 0) - pR 0 cos 0 + , 

... (4.135) 


respe ctively. 

me oonstaiits of integration are moSlfied as indicated 


= ■*■11 *■ ■‘^12 "I * ®'*13 p ’ 

^0^ 

Sl = ®ll ®®l3iy ’ 

P^a 

Ag = Agj^ pH + -^22 "I * ^ ’ 

PE 

3g = Bgj^ PE+ Bga a+ SBgS ^ ’ 
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^oR 

C, = pr’ CIR + 80^3 , 

R E^ 

Cj = Ogj^ pR® + O22 •" ®°23 ’ 

I'jjP® 

O3 = 03^ pR® + Cjz RH 8O35 . 

c, = C^, pE'^ + 0^2 ^ ^*^43 “;^’ • 


4 - -41 -t-'- * 3-; 

l-,e o=nrt.^ts of xntegratxon are determxpcd from 
the oon:x-uona g«en P7 B^s. (4 .12) trough (4.l4). 


D(x) = 0 


E(x) = 


8ar^P. 

''0 E-‘^3^' 

-rs 


r q. ^“3^ On 

■ (p «n + E - ^ 


'Sa^ 3-f 


4 0 

’2“ 


^ (Ji coo 'p — SXH cI)q) 


4E 


(L ^ 2x) + 

L 


All = ^12 = A^3 = 0 


A21 - 4a ^ag 
Aod 

^22 = V 


2af 


4 N, _ -p 

Ag3 = -(-4l"a- * 8“) "1 2ag "2 


Ir 


% IsJi - 1°-— 

“11 “ *0 *0 * 5 aJ " -laj^ag' * Sa^ 4a^ag 


1 I 


’l2 <i)r 
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Sl 


®13 " 




tan 


!3. 

2 a, 


'(f^+fg tan 4^) 


^ a I tan 

f r' 

Bgi = ^0 tan (1 + 2 i:^' - 


B 


B 


22 - J 


’23 


2^ 

0 

a„ f 

"■ ’~2a, 


3 *2- tan ♦„ + (lAr * h 


^ a ^ tan 


o 




E 


11 


°12 = °13 = ° 


'SI 


* *° ^9 

c°s *0 + ^-q 


*o( *0 - p > % 

4 Sir 


■"o ^10. 


Sa^ 


'22 


'23 


^9 ^ iVl’Jl . - sm 4, 

S-a^a-' *• 4ag SSg ° 


a. 


^ (vp)%- 4 f 

2 


~^2ar^ 4: ^ ^9 " 4:a-- ' 0 


‘•’o ^7 *’’0 % 


= 


31 


" 4 ar 


G 


G ,7 0 — 


51 

i 

0 


£l 


p, _ ^ ^7 

®33 - 4ag '■®'i®-2 9 


3 "6 


♦- fs *'0 ^6 


C 41 = ao^ *0 - 



14 1 


^42 = 


'45 


2 6 __ _ _ g ^ 

o 

1 2 


Ao 

— Cu.- 


N ^5 ‘^o - H 

fra'a; + 2.“' % - “K. + ^ 1 *0 


where f^, r = 1 , lO ore given hj 


f^ = sin Ep - 2p 


o 

(sin 2p - 2^) + :^ sm p(l - cos 2^) 


fg == COS \^i^d.j.i r^p 2 


f = sixi 2p _ 4 cos ^ cosec ^ + 2(<|. - ^) 


2 2 

f ^ = 3 cos + (cos p ^ ^ cos p) cos 2p - ■^ cos 


fg = 3 + cos 2p 

fg = 2p cos 2p + 8 cos p^(sin p - ^ cos p) - sm 2p 

= sin p - p cos p 

fg = 2 cos p 4- 2 COD 'i'o + P ■‘' “^ 0 ^ *^0 *" 

2 

f = 2 cos p ~ 2 cos <{>q- 2(4>^ ~ p) sin p-( <|)^-p) cos p 

i/ 

f = 8 (4 + S^) cos p cos p _ (5 + 2p^) cos 2p 
10 


llie stress resultonts are summarised oiter substituting 
the constants of integration in the modified form as follows : 
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(i) FOi 0 < 'S’ < 


II 


X 


(p I ■” (L-2x) ... (4.136) 

2' ~ L 


H 


X* = wk (P *0 + I - '> 

^ Jj 


... (4.137) 


N. = pPu(B cos <t> - ff* sin <l> + <t> 




11 "■■-■- • -" ■ ■ '0 ag 


■) + q. (3^2 


Fr(<}> ) 8P R a 

+ -I—) + -I- (B^3 cos ♦ - ^5 ?3(0) 


... (4.138) 


O T 

= pR (P^^ cos <|> + ^ P^(<f* ) + Cg^ _ cos 4, - (|) siTx ) 

2 r 


p (4 ) 8r 

4 - Cj_ri.(B^^ G^s 4. C22) 13 ^ 

2 -ij 


^2^4.( ^ ) 


^23^ 


... (4.139) 




Q = pR (y° r^(4' ) - 4 > cos (J) - B sin <i> ) + - S ^2 sin <}> ) 

A a. ir> ty X-i- 9 


ij 


■?- (■-XT — ®is 


» • • (4 . l4o) 


( 11 ) ror p < (^ < c|>o 


I = ^ 
X R 


SSiryJ-' 


, . 4 "'°'3-^0s 

(P *^0 ■*■ R ■• 

li 


^^"1^2 




I?5(0(x)(l-x)... (4.141) 


H 




n Sa-R. 

t 4 3 Ov 

I^P '>0 E ' 
”^pg 




J'6 


IP 

(<{>)(L-2x)+ “p"- (L-2x) 

L 

» . » (4 .142) 
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= pR(A 2 ^ sin (J) + cos <{> - <!> sin ^ ~ — - B^(4>)) 

1 2 




+ ci(Apo sin + Bgg cos (fi - “ - £■ -— ) 

1 ^ 


8p^r. 

+ — g-" 


Ag 3 son ♦ + Bgj oos ♦ + 


■) 


(4,l43) 


o ^0 

= pS'^CAg^ sin <> + cos ❖ - cos 4)- 4> sin 4> + ^ 7 ^ 


'21 


1“2 




C 3 I ^22 ^ 


8P R 

+ C 32 <> + G^g) + 


Ag 3 ^ ^ "^^23 cos y 


" ^ir;!^' + 'T^ + °33 ^ 43 , 


... (4,144) 


Q^ = pG(no^ cos - Bg^ sin 4 - 4 cos 4 + lO 


21 


° B.,-.! 4) + C^^) 


r2 


,(0 


lO ' ' < N 

+ o(A-)p cos 4 — ®92 '*’ a^a^s "** ^32 


^l-“2 


SBoR 

JJ 


3/g S'/I 

Ao^ COS <) - Bgj sm 4 - + 2 ”^" 10^‘^^**’'^33 


'■23 


1 2 


... (4.14 5) 


whe re 


r (4 ) = (2 cos ^ ^ ~ COS ^ - cos O ( COS 4 - cos ^) 
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Pg(«) 

PjC*) 

P4('>) 


Pg(*) 

I'gC'J) 

I've 

Pg(«) 


2 cos ^ ^(sxn <ii ^ 
(J) s in <> cos ^ ^ _ 2 

CCS ^ ^(2 cos <{’+<{' 

a2 cos 2<i’ 

“I cos 2^ + — 

cos p - cos <1> 


. ci> ^ sin 2<J' 
cos ^ ) + — cos 2p ‘4^-“’ 


cos 2^ cos 2’!> 


cos ^ cos ^ ^ + '““2 — 


sin <)) + <!' cos ^ cos p 


( *i>Q - <i> ) cos p - sin <j^+sin <J 


sin <) 

COD p - g — - 


<> sin (J 

^ cos p i- CQS^+ 


cos p cos <) - sin ^j) + 2 <> cos cOs ^ ^ cos 2p 

sin 

- “F" 

sin <> <) cos (j 

= (J- c os p — 4. — * — — — - 


4,7.2 Determination of limit load 

Ilie limit loads vhicli are determined using the H 
and M(^ - yield criteria are denoted by end p^g respec- 
tively, the actual lower bouna p is tne smallex of the two 
values. This is dete mined as follows : 

(i) Due to N^-yield criterion : The maximum value of in 
compression occurs at x = ^ in fhe loxioitudinal direction and 
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at <!' “ P ^ tra-iis verse 

Eq, (4,io6) il- IS oTotained 


il 


X max; 



$ 4-^ 

^0 ^ E 


c^irecxioii. Therefore from 


(cos - cos 


..» (4.146) 


The yield condition ^iven by (2.26) and. nq_, (4.i4S) give 


■Ll cos 6 


8P 




o 3 


^0 ^ 


4 . 1 

o 


... (4.147) 


(xi) lie "CO yield criterion ; As explained earlier 

in tnc Section 4.5.2 rhe failure of toe shell occurs in the 
transverse direction at ‘ere crown. Representing the values 
of IT^s and II<{, at the crovm by IT and respectively, they 


are evaluated using Eqs . (4<l38) and (4,i39). Thus 


sr.E 


a^ 


X . VUX Xt 

h’ + '^13 - ^ ’ 


... (4.148) 


xnd 


n^c = pP^^^^ll ^ ^21 - 1) + lE (-£^2 ag ^22^ 


8P R 
0 


*^^13 ^23 “ a-g ^2^» 


... (4.149) 


whe re 

cos 26 

Xi = -g + — -2 cos 6 COG 6 
^^2 = ^ ^ 1* 
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Substitution of 1 and M m the yield condition giv^n by 
Eq, (2*30) give 


^12 - 


J- [b + c§- (B 2 + ^) - ^ (B 

1 0 


M. 


12 + ^ 


8 c P^R 

+ o (B 

L 


is 

13 “ ao 


8P R' a -j 

V ^ ^23 *■ Sr ^2^J » 


M li 
o 


• • • (4 • l50 ) 


where 


D. 


R 


(B 


11 


'*'0 

^ 4 ^2 


21 


1) - c 


R 




*•0 ■" ^2 


The limit load p is the smaller of the two values 
and pj ^2 * 

4*8 Analysis of Shells 

T\/o shells with different geometric parameters are 

analysed by both elastic (26) and the proposed lov/er bound 

solutions. The ultimate load and the stress resultants are 

computed at collapse. In the limit analysis the amount of 

transverse reinforcement is taken sucn that u = 0.236. 

2 2 

a and a are equal to 150 kg/ cm and 2600 kg/cm " respectively, 
c sy 

Por estimating the collapse load for a given thickness 

of shell, ‘t%ill be tne total thickness of the shell while 

using the IT -yield criterion and the effective thickness while 
c 

using the yield criterion. 
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It 3e found from the analysis that the limit load 
increases as n decreases, 'i/'hile keeping the other parameters 
as constc-iix. Therefore n is taken as equal to 0.2 in the 
present analysis. 

4.8.1 I'Tomerical Example ± (Long Shell) : 


The geometric parameters of the shell i are ; 


= 50°, ^ =1 ^ 1 = 0.l4l7, I = 3.333, | = 100, R = 8.0 m. 

The dimensions of the edge beams are ; 2W = 2 5.0 cms, 

2H = 150.0 GDIs. The factor n is taken as 0.2. The longitudinal 

2 

steel in the edge beam is taken as 77.0 cm same a& m the 

elastic design*. Elastic analysis is done for a working load 

2 

of intensity 300 kg/m . This shell fails due to IT -yield 

L* 

criterion. The solution 2 gives the best lower bound which 

2 

IS equal to 1024.0 kg/m and hence the load facxor is 3.4l, 

The values of the reduced stress resultants at critical 
sections at collapse are tabulated in Table 4.1, The reduced 
stress resultants are plotted at critical sections for both 


the elastic and the limit analysis 
study. The variations of the stress 
shovn graphically in Eig. 4,3. 


solutions for a comparative 

y. 


resultants \nxth ^ are 

h 


*The method of design procedure is as given in the 
reference 26. 
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TaTole 4.i ; Reduced stress resultants at critical 
sections for shell l (long Shell). 
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4.8.2 lumerical Exauple 2 (Medium Shell) : 

The geometric parameters of the shell 2 are : 

♦ q = 45°, P = ■! 4>q, ^ = 2.5, = 100, R = 8.0 m. The dimensions 

of the edge "beam are i 2W = 30.0 cms and 2H = 67.50 cms . 

The factor n is tdcen as 0.2, The longitudinal steel m the 

edge heam is talcen as 58.0 cm^ same as in the elastic design. 

Elastic design is done for a working load of 4oo kg/m , The 

solution 3 gives the best lov;er hound v/hich is equal to 
2 

1920 kg/m giving a load factor of 4.8. The shell fails 
due to - Rij, criterion. The optimal value of is found 
to be 0.3381. The values of the reduced stress resultants 
at the criticsl sections are presented in Table 4.2, The 
variation of tne stress resultants with are shown graphi- 
cally in Eig, 4.4, 

Tne values of lower bounds computed from the three 
solutions for the two shells are presented in Table 4,3 
for purposes of comparison. In eacn solution collapse load p 
IS computed using both the il ^ and yield criteria. The 

smaller of these t'ro values lo the lower bound of that 
solution. In like manner the lower bounds are determined 
for the three solutions and the best lower is indicated in 


Table 4.3 
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I'alole 4,2 : Reduced streos resultants at critical 
sections for enell 2 (Medium Shell) 



R 

X 

Tq 

(x=L/2) 


«0 

~ 'V- {p 

(x=0) 

0 

-0 .420 5 

-0 .1642 

-1.1992 

0.0000 

5° 

—0 *4836 

-0.1624 

-1.1496 

0.6172 

10° 

-0.6460 

-0 *1566 

-0.9892 

1.3980 

15° 

-0.8287 

-0 . 14 5l 

-0 . 7404 

2.4321 

0 

0 

-0 .9033 

-0 .1265 

-0 .4218 

3.6635 

2 5° 

—0 .6966 

-0.1001 

-0 .090 5 

4.82 52 

CA 

0 

0 

0 .0009 

-0 .0683 

0 .192 5 

5.3809 

0 

0 

0 .8500 

-0 .0200 

0 . 1200 

4.1000 

45° 

1.4 5l8 

-0 .0044 

0 .0000 

2.5133 


fahle 4,3 : lov^^er hound solutions 


Shell 

type 

iiov/er found 
solution 1 
(p in kg/m2) 

1 lov/ef hound 
solution 2 
j (p in kg/m^) 

lower hound 
solution 3 
(p in kg/m^) 

Best 
lower 
bound 
(p in 

kg/m^) 

"i T 

c 

criterion 

criterion 

. —J 

I" 

c 

criterion 

i, 

crite- 

rion 

h ] 

c 

'crite- 

rion 

9 <p 
crit- 
erion 



long 896,0 1859*0 1024 ,0 2440,0 946,0 2640 *0 1024,0 

Shell 

(Shell 1) 

Medium 2 504.0 1044.0 2870.0 1557.0 208 5.0 1920.0 1920.0 

Shell 

(Shell 2) 




distribution of stress resultants for shell 2 




FIG.4.4 DISTRIBUTION OF STRESS REULTANTS FOR SHELL 
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4.9 Design of Shells 

In this section the design procedure for simply 
supported RC c^^lmdrical snell roofs with longitudinal edge 
beams is described* Iwo shells are designed one for each 
mode of collapse. Elastic design is also carried out for 
working loads. 

Tne design procedure is indicated as follows : 

(a) for a given T/orking load, the design ultimate load will 
be estimated using suitable load factors. Dor simplicity of 
calculations, distribution of live load can be taken as gravity 
loading . 

(b) The amount of longitudinal steel in the edge beam will 
be preassigned. Elastic design may be taken as a guide. The 
factor n can be taken as equal to 0.2. 

( c) The value of v/hich gives optimal lower bound will 
be determined in the lower bound solution 3. 

(d) Collapse load will be computed from the three lower 
bound solutions. Out of these three values the manimum value 
will be picked up and the lov/er bound solution that yields 
the maximum value of tne collapse load will be noted. 

(e) Depending on the relative values of the design and 
computed values of the limit load, a and t will be revised. 

(f) The stress resultants will be computed using the 
corresponding lower bound solution. 
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(g) 'Ihe transverse steel will be computed from Eq. (2.29) 
so as to folio ' tlie force distribution at different points 
on the shell, 

(h) longitudinal and diagonal reinforcement will he computed 
form tne distributions of U and F . respectively. 

X X Y 

In all these calculations, the stress in the steel 

will he the yield stress cr , 

sy 

To compare with the elastic design, two shells vVith 
different geometric parameters ene designed. The parameters 
of the two shells are j 

Shell 1 (Long Shell) : = 30°, ^ = | <!>^, | = 3.333 and 

R = 8.0 m. Tue edge beam dimensions are : 2\{ = 2 5.0 cms and 
2H = 150.0 cms. 

Shell 2 (Hediimi Shell) i = 45°, S = !■ ^ = 2.5 and 

R = 8.0 m. Tue edge beam dimensions are : 27 = 30.0 cms and 
2H = 67.50 cms. 

In the elastic design the tvi/o shells are designed 

2 

for a working load of 4oo kg/m • Limit design is carried 
out for an average load factor of 2.5 for hoth the shells. 
Shell 1 (Long shell) fails due to the N -yield criterion 
and the solution 2 gives the best lower hound. Shell 2 
(Medium shell) fails due to the yiela criterion and 

the solution 3 gives the best lov^er hound. 



■between the elastic and limit designs 
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Por the two shells under consideration, thiclcness 
of the shell, compressive strength of the concrete and the 
quantity of steel at critical sections are presented in 
fable 4,4, The quantity of the longitudinal steel shown 
indicates 'che steel required to resist the total longitudinal 
tension developed at the centre of the span. The transverse 
steel IS estimated at crown to resist and per unit 
length. Diagonal steel per unit length is computed to 
resist the diagonal tension equal to developed at the 

critical section which is at the support. 

Prom the foregoing analysis it can be concluded that 

the long shells fail due to the 1 -yield criterion and tne 

c 

medium snell fails due to the yield criterion. 



CHAPTER 5 


IWIIEEICAL AiTALISIS OP CXLIRREICAL SHE'LL HOOPS 
5,1 General 

In the last chapters closed fort: solutions are 
developed for the limit analysis of cylindrical shell roofs 
\;ith or without edge teams. In all the solutionsj the 
distribution of one of the stress resultants is preassigned. 
Hence the solutions are not uniq_ue • It is al?\rays possible 
to improve the lower bound by resorting to some other 
distribution. The success of the method depends mainly 
on the shill of the analyst m selecting physically appropriate 
function to represent the stress field. The closer the 
guess function is to the actual solution, presumably the 
closer IS the lower bound to the actual collapse load. 

With this In view a numerical method of analysis is 
developed. In this metnody the continuous stress fields of 
the strucxure are replaced hj^ a finite number of parameters 
by rewriting the basic differential equations of equilibrium 
in finite difference form. A linear or linearized nonlinear 
yield condition v/ritten in terms of the stress parameters is 
introduced. It is well known that any problem of plastic 
stress analysts of a structure can be recast as a problem 
in linear programming. The load parameter is maximized an 
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sucli a v/ay tliat iiov/here in the domain of the shell, the 
yieia condition is violated. This method is v^-idely applied 
to the analysis of frames and plates. Tne section that 
follo^^s gives a brief review of the literature. 

Selective literature Survey 

The variational principles of Greenberg and Prager (27) 
for the limit anal3rsis of elastic-perfectly plastic structures 
rerj_uire the optimization of a linear functional subject to linear 
constraints, Dorn (28) solved the plastic limit analysis of 
plane p m- 30 A trusses using the linear programming techni- 
q_ue. ICoopman wd Lance (29) first applied linear programming 
technique to continuous structures like pla'ces with various 
boundary conditions* Hodge and Belytschko (3o) applied non- 
linear programming technique for the limn analysis of plates. 
Botn lower and upper bounds are found out for plates of 
various boundary conditions, Grierson and Graham (3l) used 
linear programming, to compute the collapse load of plate 
structures by kinematic method. Application of linear 
programming to determine tne collapse loads for plane frames 
IS a.lGO found in the literature (o2) . 

In the literature almost all the attention is paid 
to the mathematical programming techniques applied to the 
frsuiies and plates. From the foregoing it is clear that little 
attention is paid to the analysis of cylindrical shells by 
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matlie-ia-tical programminb techniques. The mathematical progra- 
muiinc techniquesare extensively covered in the literature 
(30, 34, 35). Por completeness the method of linear programm- 
ing 13 comprehensively presented in the following section, 

5.3 linear Programming 

The metaods applied to find optiniun value of a function 
are laaoVvn as optimization or mathematical programming techni- 
ques. These methods have wide application in many areas 
including the field of engineering. Linear programming is a 
"branch of mathematical programming techniques, where in the 
oh^ective function and the constraints are linear functions 
of tue varia'oles. The brief literature review presented m 
the Section 5.2 suggests tue application of this method in the 
context of the limit analysis of structures. 

A linear programming problem can be stated as follows : 
Defining a function usually Tmovm as an objective function 
&s 

r 

2 = 2 lo X , ... (5.l) 

0=1 

it IS to be found the value of x , o = l,...,r such that 
IS maximized or minimized, x uas to sa-tnsfy the constraints 
of tue form ; 

r < 

2 X. { ^ } b , 1 = 1, 

D=1 ^ i 


j n ^ 


... (5.2) 



where only one si^n holds for each constranat. The k. , a , 
anc! h^ are assuried -to be constants. In addxtionj the value 
of chc vanxabie n should be positi-^/e. Thus 

J 

> 0 for all (5*3) 

In linear progreciuiin^^, constraints mean the restrictions 
of the form (5.2), 

5,3,1 Unrestricted Variables 


As already said that the variable x should be non- 

neyaiive. hut in most of the linear prograimung problems 

sorxie or all of the variables can take any sign. In such 

situations a variable x is expressed as difference of tv/o 

J 


varis-bles as 


^2 ^3 “ ’ 

xp >0. 


... (5.4) 


Deoending on the relative values of x' and x" , x can be 

J d J 

negative, zero or positive. But in the solution to the 

Imeor programming problem eirher x’ or x’* \7ill be present 

J 3 

(3 5) but not both simultaneously. If x' is present x” = 0 

J 0 

and vice versa. 


5.4 formulation of the Problem 
5.4,1 Finite Difference Method 

Finite difference method is used to replace the 
diCferential equa'^ions of equilibrium (2,i5) through (2.18) 
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"by al;^eToraic equations. In doing so, a second order curve 
IS assumed, to find out first and second derivatives of a 
function. If y , y^.. and y represent tnree consecutive 
fuiiction values at an interval of Al, tnen tne first and 
second derivatives a.t tne point; i are i^iven ty (24) 


dy 

dx" 


1 


SAx 


2 

^ fi-i 

dx'^ ( ^ x) 


2y + y 
•^1 -^ 1+1 


. .. (5.5) 


... (5.6) 


In case if (i+i) is a point on tne boundary, then rhe first 
derivative at tne point ( i+l) is given by 





2 Ax 


... (5.7) 


The espressions (5.5) through (5.7) are used for first 
and second derivatives in tnis thesis. 

5.4*2 Lquilibrium liquet loxis 

To generalize the netiiod. the equilibrium equations 
a^re transformed into nond miens lonal form, I'itli tms in view, 
the following transformations are made to convert the variables 
iiito nond mens lonai quantities : 



l64 


X 

L » 


R T R *^0 

r - - , t _ :|r , Sq - i,j > 


I 


n 


X 


X 


n 




iM 


X(j) 


^ = ir- 

O " 0 


■t\T 

I'i 


and n. = 


■<}> 




(5.8) 


1 ,/liere x is nondimensional ccx)rdinafce in I. direction, r and ^ 
one nondimensional shell parameters, s is the ratio of S 
and II Q, n^ and n^ are nondimensional inplrne direct forces, 
n,^^ IS the nondimensional inplane shear force and is 
nondimensional bending moment, Rurtner, if p is the 

^ ' V/ 

intensity of tne v/'orlomg load and X is a load factor then p^, 
v/liicli IS the collapse load in nondimensiona,! form, is given 


by 



* , , ( 5.9) 


Using the foregoing transformations, "clie e equilibrium 
eqao Irons (2.i5) through (2.1?) are simplified as 



= 0 , 


mm* ( 5 « 10 ) 

^ ;r- + =0 TV 

o:c 


6 il . S i'i . 

+ So r 

^ ° 5x 

4t 

6m p _ _ 

= — X — “ t s m <j> 5 

6 4 <^c 

m m m ( 5 • ll) 

, 2 

° “(J) 

n^ + — —9^ = . 

^ 4t cT 

°0 

t cos 4 , 

. . • ( 5. l2) 


Since IS not a generalized force, it is eliminated 
bj^ using the equilibrium equatxon (2,i8), 
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Partiier, tlie differential equations of equilibrium 
(5.10) through (5.12) are replaced by algebraic equations 
using tbe finite difference transloruations (5.5) and (5.6). 
Since botii the geometry -nd loading of tbe shell are symmetric, 
only a quadrant of the shell is considered for the analysis 
as shown in S'!! . 5. la. Referrinf to the Jig. 5.1b, let 
i and D locate the position of a grid point and let Ax 
anci ^ <t> be the mesh sizes m the x and f}* directions 
respectively, then Eqs. (5.10) through (5.12) at any grid point 
(1,3) are expressed as 

... (5.13) 


So 1 


8t^4> 


^<}>(i,3+l) <t>(ijD~l) 


p ^ _ 

_ X t sin i{> , 

% 


... (5.14) 


+ C%( 1,3+1) " ^“<!>(i,d) ^'^(!.( 1,3-1) 


'Hi53) 4t(A(i)) 


X t COS (j) 


• • ( 5* 15) 
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5.“- .3 Syrametry and Boundary Conditions 

Tne stress resultants n^, n^ and ui^ are symmetric 
and n.^^ is antisymmetric. Bor shells v^ith free longitudinal 
edj,es the following boundary conditions are used : 


(x) at X = 0, Lj = 0 , 

(li) at = 0. 


• . • C 5 • l6 ) 


5.4.3 I'iatrix Uotation 


At each grid point, the type of equations (5.l3) 
througii (5.15) have to be generated. As such there v/ill be 
a number of equations for the part of the shell under consi- 
deration. Bor convenience these equations are expressed in 
the 'lo-trix form. Tnis form is beneficial for computer 
ane.lysis . 


If tne mlcnown stress resultants ate cue 

8.re taken as variables and are represented bjr x 

1 

and if { X } define a column vector, bhen {X} 
of the vector expressed as 


grid points 

5 the transp ose 


{X}^ = (x^, X2,...,Xj^)* ••• (5.1?) 

let j^A^ define a rectangular matrix of size m x n 
uherc m is the number of equilibrium equaticns and n is the 
number of unknotvii stress resultants, let iS) define a 
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t ) IS a vector of size m x 1 containing tlie variables 
corresponding to tlie independent columns, 

{ } IS a vector of size (n-m) x i containing the 

varieloles corresponding to the dependent columns. 

Iluliiplying Iq. (5.19) hy on both sides 

[33-1 Cb3 Cb]-" [o 3 {xp=x Cb3-' <i>>. 

... ( 5 .EO ) 
or 

{X^} = X {E}>. j ••• (5.21) 

where 

{E> = Cb3-1{D), 

Cr3= CB3-1 [IcD. ...(5.22) 

fhus the variables in the vector are expressed 

JL. 

in terms of the vain able s in the vector 

o 

5.4: .4 Objective function an(? Constraints 

Any linear programming problem lio-s bo optimize a 
function while satisfying a set of linear inequalities or 
equalities (the constraints) of the form (5.2) and non- 
negativity restrictions of the form ( 5.3) . In the present 
problem, the load factor X is taken as an obDective function. 
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Thus 


Z 


o 


X , 


... (5.23) 


in vfliich X IS taken as a variable and its value will be found 
such uliat Z^ is optiiuun. 

The constraints are formulated such that the stress 
resultants or a combination of stress resultants should not 
exceed a specified limit at each 2,rid point. The stress 
result.ints acting at each grid point are 

Q^. Her© IS not a generalized stress and hence no restric- 
tioxi is necessary on it. 

In Chapter 2, the and yielC criteria are 

proposed. In tiie H -yield criterion, m compression zone 

IS limited to IT and tnere is no restriction on h m the 
c X 

tension zone. In ail toe closed form solutions, the magnitude 
of IT,, in the tension zone is limited by specifying the position 
of one neutral axis and satisfying the over all eg_uilibrium 
of the shell. In the presSxit numerical axialysis if no such 
restriction is made on IT , theii bhe solution becomes unbounded. 
Similstrly should also be restricted in magnitude. 

and !!(;, are restricted by bhe yield condition given by 
Eq., (2.30). These conditions give constraints at each grid 
point . 

Tnere fore the constraints at each grid point are 


obtained as 
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K. < < a- n 


X 


-I ^ 

1 c ^ 


^ ^ ^^'2 ®o’ 

£ p'* “ 5f” £ » 

*0 0 


. .. (5.24) 


vniere t-ie values ol the factors a^ and a 2 preassigned 
from pre.ctical consideration. 


The constraints (5.24) are made nondimensional hy 
dividing the first constraint by and the second constraint 
by S^, Using the relations (5.8) to express the stress result- 
ants in dimensionless form and since 

constraints (5.24) are expressed in nondimensional form as 


- 1 < . 

“ 12 ^x<J> “ ^2 ’ *** 

_b<Q(^- cn,j,<b. 


Purtner the constro^ints are modified to fit into the 
linoor programumd, problem as 

^x «i ’ 

^ “2 » 

- c n^ < b, 

- < 1 j 

Ui. 

" ^x<!> ~ “2» 

- (m^j, - c n^) < b. 


. . (5.26) 
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This t 3 rpe of constraints have been developed at 

eoch grid point. Y/nile developing the constraints, the 

metrix relation (5.2i) is used to express the variables 

(stress resultants) that are contained in the vector 

in terms of the variables contained in the vector { Z }. 

s 

Thus all the constraints are expressed in terms of the 
variables in the vector {Z_}. This clearly shows that the 
number of variables and consequently, the size of the matrix 
IS reduced. If no such elimination is done as given by the 
matrix relation (5.2i), the equilibrium equation given by 
tne na.trix relation (5.l8) can be treated as equality cons- 
traints, which requires additional computer memory. \/hile 
treating equality constraints, artificial variables are to 
be added to form tiie initial basis, Further, since the load 
factor X IS taken as a variable, \vhen the vecior { D } is 
transferred to the left side, the rignt hciid side contains 
zero vector in the relation (5.l8). If bhe equality constraints 
contain zeros on t-ie right hand side, artificiel variables 
enter the basis at zero level without improving the objec- 
tive function Yhiich becomes a degenerate problem. This 
process will be continued till all the artificial variables 
are eliminated from the basis. Hence convergence becomes 
poor because of the roundoff errors. 

After formulating the constraints, they can be 
expressed in the matrix form as 
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C&3 ^ Xg} < {Y}, ... (5.27) 

where the matri:: ^ C [] is formed by the coefficients of the 
constraints ond tbe vector { Y } represents the right hand 
side of the constraints given by inequalities (5.26). 

Ihus the objective function and the constraints are 
given by the relations (5.23) and (5.27) respectively for 
the problem of linear programming. 

5.0 Computer Programme and Results 

A compucer programme has been developed in which the 
following aspects are dealt with, viz., (i) generation of 
the elements in tne various matrices Lb:, Lc], 

and vectors {P} and {B}, ( ii) identification of 
the variables in tlie vectors } and { X }, (m) generation 

X S 

of tne elements in the ciatrin ^ and toe -s/ector {Y }. 

Ihe revised sim'>lez method (35) is used to solve the 
linear programming problem. Tiie constraint matrix Cg-]] and 
the vec'cor{Y} in the matrix relation (5.27) are given as a 
dabs to tne revised simplex method computer programme. Ihe 
inequalities in the relation (5.27) will be converted into 
equalities after adding slad: variables. The values of the 
obi active function tne variables in the vector {X } and 
xne slack; variables are given as an output in tne computer 
program lie . 
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The values of tne variables in the vector { X ) 

r 

ar® deterrmnea from tne matrix relation (5.2i) by hack 
suhst itution , 

Thus the optimal value of the load factor X and the 
stress resultanto at different grid points a.re evaluated. 

5.5.1 N ume r i ca,! Zxa,_iD le 

A cylindrical shell 'yi'th free longitudinal edges is 

analysed to illustcate the numerical technique. The geometric 

parameters of bhc shell are ; <i>^ = 45°, R = 8.0 m, g = 2,5 

and t = 8,0 cms , In tue analysis and in the equatioxi 5,9 

2 2 

are xaken as equal to 300 kg/m and I50kg/cm respectively. 

To apply the finite difference scheme, the shell is discrebised 

4 > 

as a grid v/it:h me or size of A 5c = and A41 = as shovm in 
Pig. 5.1a. The facbors a-, And ao am "tue constraints (5.26) 
are taken as 1.5 end 3*0 respectively iii tne present problem. 

In Section 3.8.2 a shell v/ith tne same geometric 

parameters is anolysed. Tne loner hound solution 3 is bhe 

best solution and the load factor 13 equal to 1.0 50 . The 

present numerical solution given a load factor of 1.9 72. The 

reduced stress resultants that are obtained a.! different grid 

points by the nuiierical solution are given in Table 5,1. 

The distribution of stress resultants by both the analytical 

and numerical mebnods is shown in Rig. 5.2 for purposes of 

2 

comparison. distribution is drawn for x = g- L and for 
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Table 5.i ; Reduces stress resultants at 
different grid points 


.d point 
Ro. 

IT 

-v* 

So 

'•o 

*'6 

1 

-0 .2545 

0.0- 

-0.3472 

-1.6400 

2 

-0 .50 

0.0 

-0 .4303 

-1.7444 

3 

o 

• 

o 

0.0 

-0.3108 

-1.7724 

4 

0.0 

0.0 

0 

• 

0 

0.0 

5 

to 

• 

o 

1 

0.0 

-0 .33 50 

-1.6296 

6 

-0 .3140 

-0.8100 

-0.3377 

-1.6290 

7 

0.0 

-1.6168 

-0 .3230 

-1.6288 

8 

0.0 

-2.3435 

0.0 

0.0 

9 

0 .1275 

0.0 

-0.1293 

-0.6900 

10 

0 .2 500 

-0 .9875 

-0.1378 

-0.7036 

11 

0.0 

-1.9738 

-0 . 1400 

-0 .7200 

12 

0.0 

-3 .0 

0.0 

0 

• 

0 

13 

0.2120 

0.0 

-0 .0600 

-0 .2 550 

14 

0 A ‘±00 

-0.89 50 

-0 .0640 

-0 .2480 

15 

0.0 

-1.415D 

-0 .0653 

-0 .24 50 

l6 

O 

• 

o 

-2.70 

0.0 

0.0 

17 

0 .497 5 

0.0 

0.0 

0.0 

18 

1.00 

0.0 

0.0 

0.0 

19 

0 .0 

0.0 

0.0 

0.0 

20 

0.0 

0.0 

0.0 

0.0 



“Lower bound 



FIG. 5.2 DISTRIBUTION OF STRESS RESULTANT 
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it IS drawn at support. Iv^ and distributions are drawn 
for X = ^ . 

5.6 Discussion of Results 

Hurierical analysis is carried out for the shell with 

two different mesh sizes. In the first one ax = ^ and 

4>o _ ^ 

A(j> ~ and m tne second mesh ax = ^ and A(fi = ^ are 

considered. Values of i.2io and 1.972 for x are obtained 

in the first and second case respectively. 

In the present numerical analysis double precision 
accuracy is necessary in computer programme. It is observed 
that the single precision value for the load factor is about 
1.5 times the value obtained with double orecisioix accuracy. 
Because of the inadequacy of computer storage, further division 
of the shell into finer mesh is not carried out to test for 
the convergence of this technique. 

Rrom the values of the stress resultants presented in 
Table 5.1, it is observed that tne longitucl inal variation of 
M. and N, at a oarticular angle 4> is not considerable. This 

ifi " 

indicates the sane trend as obbained in tne closed form 
solutions wherein it is found that and are independent 
of X. rurther the maximum value of obtained in the present 
numerical solution is about i,75 tries that of the value 
obtained by the lower bound solution 3 (Pig. 5.2c). 
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The distribution of (Pig* 5.2a) is linear m hotli 

-X. 

the compression and tension zones v/hereas in the lower hon-n d 

solution 3 it IS parabolic in the compression zone. In addi- 

2 

tion, maximum values for 1 occur at x = ^ I in the present 

X o 

numerical solution instead of at centre of span in solution 3. 

And also it is found N becomes zero at all the grid points 

5 

lying on the transverse cross section at x = ■g“ I< * 

Pig. 5.2b indicates the distribution of as per 
the numerical solution and lower bound solution 3. 

Prom Pit- . 5. 2d it can clearly be found that there is a 
wide discrepancj^ in tne magnitudes of 11(1) by both numerical 
and loiier bound solutions. The value of !!({> obtained by 
numerical techni'‘ue is quite high and it is 7.6 times that 
of the value obba.inod by lower bound solution 3 at crown. 

Prom the roregoing discussion it can be concluded that 
the stress resulcrnt IT^ is quite sensitive in the present 
numerical solution. The solutions are converging as the mesh 


size decrea^ses. 
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ae oiEwim 

6.1 Suxau-ar^ and Discussion 

Cnis tiiesis is mainly concerned w ith the development 
of lower hound solutions for EG circular cylindrical shell 
roofs. Dor the shells under consideration, the ends are 
simply supported and the longitudinal edges can he either 
free or v/itli edge beams. The analysis is carried out for 
uniform gravity loading. 

limit analysis requires a yield condition \;hich 
characterises the failure of a material of the structure. 

To the best of the author’s knowledge no yield condition is 
available in the literature for reinforced concrete cylindrical 
shells. Doca-use of this fact, tv/o independent yield criteria 
are proposed viz., ( i) E^-yield criterion, (ii) yield 

criterion. Tne IT^-yield cricerian states that a snell fails 
if reaches a specified value m compression. The 
yield criterion states that a shell fails due to the inter- 
action of bhe transverse moment and membrane direct force 
. This yield condition is eq_uivalent to the failure of a 
reinforced concrete section subjected to a moment II and an 
axial force IT (22) . To reduce the complexity of the problem 
the yield curve of the yield criterion is suitably 



Ixnearized by inscribing a polygon# lliis modification is on 
the consej-Tot I’/e ='ibe. A shell fails if at least one of the 
two yield conditions is satisfied. 

In the r -yield criterion the failure of a shell is 

initiated, due to crushing of concrete in compression, This 

may cause a sudden collapse of a shell. To avoid the sudden 

collapse, the limiting strength of a shell in compression (W ) 

c 

can be reduced by multiplying with some preassigned factor. 
Because of this reduction of iJ , a reduced collapse load is 
obtained# But in an actual collapse the shell fails due to the 
yielding of steel in the tension zone before it actually fails 
due to ccishing of concrete. 

In the present lower hound solutions, H is not 
restricted in magnitude# Diagonal steel vjill be provided to 
resist diagonal tension developed due toH^^# In some 
occasions \'here the computed diagonal steel is heavy, a 
restriction can be placed on the magnitude of • G5pde 
specifications regarding the spacing of the reinforcing bars 
will flic the magnitude of il . in such situations. 

X 9 

Tiiie drawing the yield curve forK^-N^ yield 
criterion, the value for reduced reinforcement ratio 

A Q 

— JJ) xs tsdeen as equal to 0.236 in Eq. (2#29). This 
quantity of si eel is the value required for a balanced 
rectangular section under pure flexure as per 1 .3 . Code (25) 
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lay ultifiic.te strength design. The contents h and c vfhich define 
the straigat line equation given by Eq. (2.30) are computed 
as 0.833 and 2.25 respectively. For any other value of y, 
the yield curAre can be dravm using Eq. (2,29) and the corres- 
ponding values of b and c in Eq. (2.30) can be determined. 

Three lower bound solutions are developed for shells 
with and \7ithout longitudinal edge beams. In all tnese 
solutions, the distribution of is preassigned in the 
transverse direction and the other stress resultants are 
solved using the differencial equations of equilibrium, 
boundary and continuity conditions. 

Tor shells with edge beams the transverse moment 
at the junction of the shell and edge beam is tahon os zero. 

And also it is considered tnat the edge beam cannot lesist 
horizonto-l forces. These assumptions are the same s,s m the 
elastic anal;'’sis (37). The distribution of lon_, itulinal force 
across a cross section of the edge bean is tadoen as trapezoidal 
(fully in tension). To determine the collapse loaZ the edge 
beam reinforcement has been preassignod. Elastic a.esign can 
be talcen as a guide. The factor n specifies the distribution 
of force in the edge beam (Fig. 4.ic) . It is observed that 
as n decreases, collapse load increases. A minimum value oa. 

0.2 for n is taken m the present study. 
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In all thG lov/er bound solutions the angle ^ which 

defines tno position of neutral axis is preassigned. It is 

found fron bne analysis of shells as the angle p, increases 

collapse load increases. Therefore the value of the angle ^ 

be limited such that the longitudinal steel can be accommodated 

in the tension zone <J>q) * In the present study ^ is 

2 

taken as eaual to ^ <J).. 

o o 

In the lower bound solution 3 (Sections 3.6 and 4.7) 
the angle r virhich gives the position of maximum value of 1 
in the compression zone (O < 7 < p) is tai;en as a variable 
while other parameters of a- shell are treated as constants, 

For shells ’'itli free longitudinal edges the optimal values of 
which give the best lower bounds are presented in Table 3.1. 
For shells with longitudinal edge beams the optimal value of 
depends on the factors namely depth of the edge beam and 
the amount of longitudinal steel in the edge beam. Therefore 
the optimal value of p ^ has to be determined for each problem. 

For shells with free longitudinal edges the variation 
of nondimensional load p^ \/ith jj' ratio for different values 
of are presented by way graphs (Fig. 3.2), Tnere graphs 
are useful for preliminary design of cylindrical shell roofs 
without edge beams. From the graphs (Fig. 3.2) it can be 

•p 

observed tiiat for ^ ratio greater than some value j shells fail 
due to bliG yield criterion and the collapse load is 
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R R 

constant irrespective of ^ ratio, fhis limiting “ ratio 
beyond 'vfricli collapse load is constant v^ill have different 
values for shells v/ith different semicentral angles (4*^). 

Tno distribution of stress resultants for a long 

shell vfitu free longitudinal edges is presented in Rig. 3,3. 

Ihis shell nas failed due to the N -yield criterion. And 

c 

also It is shov«i on Rig. 3.4 the distribution of stress 

resultants Tor a medium shell which has failed due to the 

-Ijj, yield criterion. Similarly the distribution of stress 

resultants for long and medium shells with longitiidmal edge 

beams is sho\ai in Rig. 4,3 and Rig. 4,4 respectively. Railure 

of the long snell is governed by the R -yield criterion and 

c 

that of the medium shell is governed by the yield 

criterion. Rurther there grapns give an idea about the 
distribution of forces am collapse of a. shell, 

Ror comparison, t\/o shells with free lor.t^ itudinal 
edges (one long and one medium) are designed by both the 
present limit and elastic analyses. The results are presented 
in Table 3,5, It can be observed that there is a considerable 
economy due to reduction of concrete strength, tniclmess of 
shell and also quantity of steel as compared to the elastic 
design. Besides a comparative design statement for shells 
V 7 ith longitudinal edge beams is also presented in Table 4.4. 
For these shells therb is no saving of materials but in fact 
there is an increase of reinforcement as compared to the 
elastic design* 
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111 Clicpter 5 an attempt is made to introduce a 
numsri-cml tecnnique to obtain lower bound solution for a 
EC cylindrical shell roof, A comparison is made betvreen the 
distribution of stress resultants as obtained by both numerical 
and closed form solutions. 

6,2 0 on clu s ions 

Ivfo independent yield criteria ■viz.,(i) IT —yield 

c 

criterion, (ii) yield criterion are proposed for EC 

cylindrical shells (Section 2.3.6). These approximate yield 
criteria can be adopted toll an accurate yield criterion is 
developed supported by experimental evidence. However the 
lower bound solutions that are developed in the Ciiapters 3 and 
4 are valid even for a revised yield criterion. 

The lower bounds computed from the proposed limit 
analysis coincides with the upper bound evaluated, for a beam 
mechanism, talcing the shell a beau of tnin curved cross 
section, If the shell faj-ls duo to f -yield criterion (See 
appendix) . hence there solutions give the enact limit load 
If the shell fails due to the H -yield criterion. 

Prom the analysis it is observed that long shells 

fall due to the E -yield criterion. Por medium shells v/ith 

c 

higher values of semicentral angle ( (fi^) say 45° and above the 
yield criterion governs the collapse. 
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Lo\/Gr iDouaad solution 2 always gives a better value 
than the solution i. Lower bound solution i may be used for 
preliminary design for the simple reason that the expressions 
for the stress resultants are fairly simple as compared to the 
other t '■Q solutions. 


In general, it can be said from the Fig. 3.2, solution 2 

gives tile best lower bound if a shell fails due to the H -yield 

c 

criterion. Further, solution 3 gives the best lower bound if 
a shell fails due to the yield criterion. 


The stress resultants 11 , Q,, and 11 . are neglected 
m the preseiit study. This assumption is satisfactory if the 
span L exceeds its chord width B by i.Q or over (24), Hence 
these solu'cions are satisfactory for shells satisfying the 
ab 0 V3 1 im it c t ion . 


Tne numerical technique developed in Chapter 5 serves 
as an alternative to the analysis of a SO cylindrical shell 
roof at coj.laose. It avoids the difficulty of selecting 
various statically admissible stress fields for complicated 
structures lil^e shell roofs. The stress resultant is very 
sensitive in tne present numerical solution. 
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6*3 Pur-tiier Scope 

It 13 proposed : 

Cl) to develop an accurate yield criterion for EC 

cylindrical shell with eyper imental evidence, 

(2) to develop a letter lower "bound solution for 
c 3 "lindrical shell roofs, 

(3) to analyse short span cylindrical shell roofs, 

(4) to analyse multibarrel shell roofs and 

(5) to compare the numerical results by the finite 
element approach. 
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Upper Bound Solution 

Bor completeness, the upper hound solution lor simply 
supported RC circular cylindrical shell rool evailahle in 
the literature (22) is presented her®. 

let IT* denote the force/unit length in tension, zone 
P 

of a shell. The shell is considered as a heam of curved cross 
section. The load per unit length P of the shell under a 
uniform grarity loading of intensity p is given as 

P = 2pR ’{’q • •••(!) 

Assuming a simple beam mechanism v/ith a central 
plastic hinge, the external T/ork W-g. is 



= pREi a j « * , (2) 

where A is an arbitrary deflection of the shell at its 
centre during collapse. 

The plastic moment capacity M , of the shell as a 

y 

beam with uniform compression at collapse is 

= 2U^ R((^q ~ a, ... (3) 

where a is the Xever arm (vertical distance between the 
centroids of the compression and tension) and is given as 
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R( sin ^ - pi sin <i>^) 


... ( 4 ) 


Tlie internal work Wj, during the plastic collapse 


HJ A 


Y/x = 


I - I 


• • * (5) 


Using the relations (3), (4) and (5) 


= 


81Y’ R A ( ({)|^ sin p - p sin 4'^) 


. ...( 6 ) 


Equating tne external work done YYp and the internal work 


Uj It IS oh tamed 


8N' R(‘i>Q smp - p sm <{) ) 


t 




where is the upper hound. 

Prou the static equilihrium condition of the shell 
( algebraic sum of the forces in the longitudinal direction is 
zero) it rs found 


p TtT-TT ’ 


... (s) 


where 1 = a t. 

c c 


Substitution of Eq. (8) m Eq. (7) yields 
8R^R ((j)^ sm p - p sm ‘J’^) 


m * 0 (9) 



wiiicli is 3 c" -le as the lo'v/er hound given hy Eq. (3.42), 


Siailarly it can he shown for shells with longitudinal 
edge beam also. Further it can he shown that the lov/er and 
upp-er hounds are same foi other lower hound solutions also 
when the shell fails due to the H -yield criterion, 

Trora the foregoing it is concluded that an exact 
value for rhe collapse load is found when the U —yield 
criterion governs the collapse of a shell. 
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